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PREFACE 


The  High  Altitude  Effects  Simulation  (HAES)  Program  sponsored  by 
the  Defense  Nuclear  Agency  since  the  early  1970  time  period,  comprises 
several  groupings  of  separate,  but  interrelated  technical  activities,  e.g., 
ICECAP (Infrared  Chemistry  Experiments  - Coordinated  Auroral  Program).  Each 
of  the  latter  have  the  common  objective  of  providing  information  ascertained 
as  essential  for  the  development  and  validation  of  predictive  computer  codes 
designed  for  use  with  high  priority  DoD  radar,  communications,  and  optical 
defensive  systems. 

Since  the  inception  of  the  HAES  Program,  significant  achievements 
and  results  have  been  described  in  reports  published  by  DNA,  participating 
service  laboratories,  and  supportive  organizations.  In  order  to  provide 
greater  visibility  for  such  information  and  enhance  its  timely  applications, 
significant  reports  published  since  early  calendar  1974  shall  be  identified 
with  an  assigned  HAES  serial  number  and  the  appropriate  activity  acronyn 
(e.g.,  ICECAP)  as  part  of  the  report  title.  A complete  and  current  bibli- 
ography of  all  HAES  reports  issued  prior  to  and  subsequent  to  HAES  Report 
No.  1,  dated  5 February  1974  entitled,  "Rocket  Launch  of  an  SWIR  Spectrometer 
into  an  Aurora  (ICECAP  72),"  AFCRL  Environmental  Research  Paper  No.  466,  is 
maintained  and  available  on  request  from  DASIAC,  DoD  Nuclear  Information 
and  Analysis  Center,  816  State  Street,  Santa  Barbara,  California  93102, 
Telephone  (805)  965-0551. 

This  is  the  jsixty  second  report  in  the  HAES  series  and  covers 
technical  work  performed  under  DNA  contract  001-76-C-0138  from  1 November 
1975  to  31  March  1977.  The  main  purpose  of  the  work  reported  herein  is  to 


interpret  and  analyze  data  obtained  from  the  ICECAP  program  so  that  models 
can  be  verified  and/or  constructed  to  describe  the  nuclear-disturbed 
environment . 

We  would  like  to  acknowledge  the  full  cooperation  of  numerous 
people  involved  in  the  ICECAP  and  EXCEDE  programs  and  their  aid  in  freely 
making  available  to  us  their  data  prior  to  its  publication  elsewhere. 
Special  thanks  go  to  Mr.  James  Ulwick,  Dr.  A.  T.  Stair,  Mr.  Ned  Wheeler, 

Mr.  Robert  O'Neil,  Mr.  Philip  Doyle  (all  of  AFGL) , Mr.  Richard  Hegblom 
(Boston  College),  Mr.  William  Grieder  (Visidyne,  Inc.),  Dr.  Irving  Kofsky 
(Photometries),  Dr.  Kay  Baker  (Utah  State  University),  and  Dr.  Murray  Baron 
(SRI).  The  encouragement  of  Dr.  Charles  Blank  (DNA)  is  also  appreciated. 
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SECTION  1 
INTRODUCTION 


This  is  the  last  in  a series  of  reports  describing  work  performed 
at  Mission  Research  Corporation  in  support  of  DNA's  field  programs  in  Alaska, 
especially  ICECAP  and  EXCEDE.  It  extends  the  analyses  described  before 
(References  1 through  4)  by  including  the  results  of  studies  on  two  new 
auroral  events,  and  presenting  additional  work  related  to  the  earlier  events. 

The  main  purpose  of  the  work  has  been  to  understand  and  interpret 
the  field  data  from  an  auroral ly-excited  atmosphere,  especially  data  related 
to  short  wavelength  infrared  radiation  (SWIR) , so  that  better  models  could 
be  constructed  to  describe  the  nuclear  environment.  In  this  endeavor,  we 
have  availed  ourselves  of  ICECAP  data  from  five  different  auroral  events,  and 
and  have  made  extensive  use  of  the  ARCTIC  code  that  was  devised  specifically 
for  auroral  calculations  and  comparisons  with  the  field  data. 

An  attempt  is  made  to  summarize  the  main  results  and  conclusions 
of  our  work,  especially  as  it  impacts  on  SWIR  radiation  processes,  and  to 
relate  it  to  the  nuclear  environment. 

Sections  2 and  3 describe  results  and  conclusions  from  our  study 
of  the  auroral  events  of  12  March  1975,  and  10  March  1975,  respectively. 

Section  4 describes  the  modifications  introduced  into  the  ARCTIC 
code  that  are  necessary  for  the  correct  calculation  of  the  electron 
temperature . 


Section  5 presents  the  results  of  our  efforts  to  infer  NO 
profiles  for  several  different  auroral  events. 


Section  6 contains  a summary  of  the  principal  results  of  our 
analyses  of  five  different  auroral  events. 


Section  7 gives  a brief  discussion  of  work  done  in  support  of 
DNA's  satellite  communications  effort,  especially  the  Wideband  Satellite 
experiment. 

Three  appendices  are  included.  Appendices  A and  B are  in  support 
of  Section  4,  and  give  details  related  to  the  heating  of  the  electron  gas 
by  high  energy  (>1  keV)  and  low  energy  (£2S  eV)  electrons,  respectively. 
Appendix  C describes  a proposed  method  for  measuring  the  NO  concentration 
in  an  auroral  environment. 
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AURORAL  EVENT  OF  12  MARCH  1975 
(Sergeant  Hydac  IC  519.07-lB) 


PRELIMINARY 

On  12  March  1975,  at  0748:10  UT,  a multi -instrumented  Sergeant 
Hydac  rocket  (IC  519.07-lB)  was  launched  from  Poker  Flat,  Alaska  into  an 
active  aurora,  reaching  an  apogee  of  about  184  km  at  225  sec  after  launch. 
On-board  instrumentation  that  obtained  data  included  particle  detectors 
(scintillator,  electrostatic  analyzer),  a CVF  spectrometer  (SWIR) , vertical 

o o o 

radiometers  (4.3,  5.3  ym) , a 3914A  photometer  (vertical),  3914A  and  5577A 
photometers  (horizontal),  and  a mass  spectrometer. 

Ground -support  data  for  the  event  were  also  obtained  from  All-Sky 
cameras,  meridian  scanning  photometers,  and  the  SRI  incoherent  backscatter 
radar. 


The  rocket  data,  along  with  what  ground- support  data  we  could 
assemble,  were  used  as  the  basis  for  making  an  ARCTIC  code  calculation  for 
the  event.  The  following  subsections  give  a brief  description  of  the  event, 
of  the  calculations,  and  of  the  comparisons  with  field  data. 

QUALITATIVE  DESCRIPTION  OF  EVENT 

The  region  traversed  by  the  rocket  was  part  of  a system  of  arcs 
that  originated  north  of  Poker  Flat  (Reference  5) , but  within  20  minutes 
before  the  rocket  launch  had  moved  rapidly  southward  (Reference  6) . Study 
of  the  Lockheed  scanning  photometer  data  taken  at  Chatanika  (Reference  7) , 
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the  scintillator  data,  and  the  All-Sky  camera  pictures  from  Poker  Flat, 
leads  to  the  following  qualitative  picture  of  the  event. 

At  the  time  of  rocket  launch,  an  arc  complex  was  located  approxi- 
mately 40  km  to  the  north  of  Poker  Flat.  The  structure  changed  in  width, 
appearance,  and  position  as  multiple  arcs  of  varying  intensities  merged, 
decayed,  reformed,  etc.  Figure  2-1  shows  a piece  of  the  data  inferred  from 

O 

the  Lockheed  4278A  meridian  scanning  photometer  for  the  time  interval 
spanning  the  rocket  flight.  The  photometer  scanned  from  north  (N)  through 
zenith  (Z)  to  south  (S)  and  back  through  zenith  to  north  again  in  a repe- 
titious pattern.  The  main  arc  complex  is  located  a few  degrees  to  the  north 
of  zenith  in  each  scan.  The  times  of  rocket  launch,  apogee,  and  impact  are 
indicated  in  the  figure.  The  changing  width,  structure,  and  intensity  of  the 
the  arc  complex  are  apparent  in  the  seven  times  the  structure  was  scanned 

it 

in  the  sample  shown  in  Figure  2-1. 

Figure  2-2  shows  particle  data  from  the  rocket-borne  scintillator 
(Reference  8)  that  responded  to  flux  from  electrons  with  energy  >4  keV. 

Note  the  sharp  increase  in  flux  as  the  rocket  ascended  from  100-  to  130-km 
altitude.  We  believe  this  occurred  as  the  rocket  penetrated  the  southern 
edge  of  the  arc  complex,  reaching  the  intense  inner  arc  at  about  120  sec 
after  launch.  The  smaller  flux  values  at  100  km,  however,  are  probably  due 
partly  to  degradation  of  the  incident  flux  at  this  lower  altitude.  As 
Figure  2-2  indicates,  the  flux  of  electrons  with  E>4  keV  encountered  by  the 
rocket  as  it  traversed  the  arc  complex  was  quite  irregular.  Exit  from  the 
northern  edge  appears  to  have  occurred  near  apogee  and  the  rocket  then 
descended  through  a region  of  low  particle  precipitation. 

Figure  2-3  shows  a pictorial  representation  of  the  rocket  trajectory 
with  superimposed  arc  complex,  as  inferred  from  the  description  above.  The 


* In  Figure  2-1  another  strong  arc  can  be  seen  north  of  Fort  Yukon  (70°  north 
from  zenith  at  Chatanika  (near  Poker  Flat))  at  about  0753  UT. 
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Rocket  Altitude  (km) 
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Figure  2-3.  Pictorial  representation  in  magnetic  meridian  plane  showing  rocket 
(IC  519.07-lB)  trajectory  and  approximate  location  of  arc  complex. 
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location  of  Poker  Flat  (PKF)  and  Fort  Yukon  (FYU)  are  also  shown.  Precipi- 
tation occurred,  of  course,  over  the  entire  trajectory,  but  it  was  most 
intense  in  the  region  indicated  in  Figure  2-3. 

Figure  2-4  (from  Reference  9)  shows  the  energy  deposition  rate  due 
to  electron  precipitation  and  to  Joule  heating  during  and  for  about  one 
hour  prior  to  the  rocket  flight.  The  relatively  low  deposition  rate  is 
clearly  evident  until  the  arc  system  moved  into  the  region  after  about 
0740  UT.  Predosing  of  the  rocket  airspace,  therefore,  would  not  seem  to 
have  been  very  significant. 

ARC  AND  EXO-ARC  ENVIRONMENT 

A quantitative  description  of  the  arc  system,  including  its  north- 
south  width  and  motion  is  difficult  to  obtain  from  the  scanning  photometer 
data  available  to  us  from  Chatanika  and  Fort  Yukon.  Camera  malfunction  at 
Ester  Dome  precluded  data  from  that  station.  Plots  of  peak  intensity  as 
functions  of  elevation  angle  obtained  from  the  two  stations  lead  to  arc 
positions  and  motions  that  are  not  entirely  consistent.  However,  if  one 
assumes  that  the  arc  system  was  fairly  stationary  while  the  rocket  traversed 
it,  the  scintillator  data  (Figure  2-2)  implies  that  the  system  extended  for 
approximately  50  km  in  the  magnetic  meridian  plane.  At  the  time  of  rocket 
traversal, ' its  position  is  approximately  as  indicated  in  Figure  2-3. 


The  spectral  flux  of  precipitated  electrons  with  energy  in  the 
range  from  about  1 to  30  keV  was  measured  with  the  rocket-borne  electro- 
static analyzer.  Figure  2-5  shows  several  measured  spectra  corresponding  to 
regions  inside  and  outside  the  arc  complex.  The  numbered  spectra  correspond 
to  times  and  altitudes  indicated  on  the  scintillator  plot  in  Figure  2-2.  Thus, 
spectra  (T)  and  correspond  to  the  region  just  south  of  the  arc  where  the 
flux  is  lower;  spectra  (3)  ^nd  (T)  are  typical  of  the  flux  at  the  maxima  and 
minima,  respectively,  in  the  arc  complex;  and  was  measured  at  270  sec 
on  the  downleg  of  the  trajectory  in  a region  north  of  the  arc  system. 
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Figure  2-4.  Energy  deposition  rate  during  and  prior  to  flight  of  Multi  rocket 
IC  51 9. 07-1 B,  as  inferred  from  Lockheed's  three-beam  photometer 
(Reference  9). 


ARCTIC  CODE  CALCULATIONS 
Model  Flux  and  Environment 
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As  described  in  References  2 and  4,  the  procedure  used  in  making 
previous  ARCTIC  code  calculations  has  been  to  deposit  the  arc  flux  in  a pre- 
dosed atmosphere  resulting  from  30  minutes  of  "drizzle"  bombardment.  The 
drizzle  flux  in  turn  was  assumed  to  be  deposited  in  an  initial,  high -latitude 
atmosphere  that  was  described  in  Reference  2,  p.  28. 

In  the  present  case  it  was  felt  that  we  had  insufficient  informa- 
tion to  define  adequately  the  drizzle  flux.  Consequently,  the  arc  flux 
was  deposited  directly  into  the  initial  high-latitude  atmosphere,  and  no 
attempt  was  made  to  calculate  an  exo-arc  drizzle  environment. 

Figure  2-6  shows  the  electron  flux  spectrum  input  to  the  ARCTIC  code. 
It  is  spectrum  @ from  Figure  2-5  (corresponding  to  one  measured  at  132  sec 
when  the  rocket  was  at  145-km  altitude  on  ascent)  with  an  arbitrary  extrapo- 
lation at  the  low  energy  end  to  0.1  keV.  Pitch-angle  isotropy  in  the  down- 
ward hemisphere  has  been  assumed  because  this  is  consistent  with  the  scintil- 
lator data  over  the  range  measured  from  about  33°  to  77°.  The  energy  input 

-2  -1 

rate  to  the  atmosphere  from  this  spectrum  is  about  52  ergs  cm  sec 


The  energy  was  deposited  in  a grid  of  25  cells  extending  in  the 
vertical  direction  from  89-  to  199-km  altitude.  Table  2-1  shows  the  assumed 
deposition  atmosphere,  with  species  concentrations  in  units  of  cm  ^ and 
temperature  (electronEneutral)  in  °K.  The  first  two  columns  are:  cell 
altitude  (km),  and  cell  thickness  (km),  respectively.  The  NO  profile  assumed 
is  the  high  latitude  one  shown  in  Reference  2,  Figure  2-11. 

The  energy  deposition  rate  in  the  calculational  volume  between  89 

-2  -1  -2  -1 
and  199  km  is  41  ergs  cm  sec  . Of  the  remaining  11  ergs  cm  sec  , part 

is  deposited  in  the  atmosphere  above  199  km  and  the  rest  is  backscattered 

to  space. 
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The  calculated  energy  partition  and  resulting  production  rates  of 
species  were  input  to  the  chemistry  subroutine  that  was  then  run  for  times 
to  1800  sec.  A great  deal  of  output  is  obtained  from  each  running  of  the 
code.  Included  are  volume  production  rates  of  over  200  states  of  excitation, 
species  concentrations  and  temperatures  as  functions  of  time,  and  volume 
production  rates  of  photons  in  numerous  optical  and  IR  bands.  However, 
only  a small  fraction  of  this  output  will  be  shown  since  most  of  it  cannot 
be  compared  directly  with  field  data. 

Results  and  Comparisons  with  Data 

The  calculated  energy  deposition  rate  profile  is  shown  in  Figure  2-7. 
Due  to  the  softness  of  the  energy  spectrum  (see  Figure  2-6) , the  peak  of  the 
deposition  is  seen  to  occur  near  the  relatively  high  altitude  of  120  km. 

The  inflection  point  at  100-km  altitude  arises  from  the  secondary  peak  in 
the  electron  flux  spectrum  at  20  keV. 

Figure  2-8  shows  selected  volume  emission  rates  calculated  for  the 
arc.  The  results  for  0,  NO,  and  NO^  are  applicable  at  both  100  sec 
and  1000  sec.  Quenching  of  the  4.3-ym  emission  from  NO^  is  not  included. 

ELECTRON  DENSITY 

Unfortunately,  no  electron  density  data  from  the  rocket  are 
available  since  both  the  plasma  frequency  probe  and  the  retarding  potential 
analyzer  failed.  The  only  direct  ionization  data  we  have  are  from  the  SRI 
incoherent  backscatter  radar  that  was  in  a meridian  scan  mode  during  the 
rocket  flight.  From  the  resulting  electron  density  contour  maps  (Reference 
6)  we  have  inferred  the  ascent  and  descent  profiles  shown  in  Figure  2-9, 
Plotted,  for  comparison,  is  our  calculated  profile  based  on  the  arc  spectrum 
shown  in  Figure  2-6.  Agreement  between  the  calculated  and  measured  (ascent) 
profiles  is  seen  to  be  good  up  to  120  km,  but  the  two  curves  diverge  above 
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Figure  2-7,  Calculated  energy  deposition  rate  in  auroral  arc 
(Sergeant  Hydac  IC  519.07-lB). 


\li 


1.  Nj  {3914A)  \\ 

2.  0 (5577A)  >\ 

3.  0 (6300A)  \ 

4.  CO2  (4.3  pm-100  sec) 

5.  CO2  (4.3  ym-IOOO  sec) 

6.  NO  (5.4  ym) 

7.  NO  (2.7  ym) 

8.  NO^  (4.3  ym-unquenched) 


no  130 


150  170 


Altitude  (km) 


Figure  2-8.  Calculated  volume  emission  rates  in  auroral  arc  of  12  March  1975 
(Sergeant  Hydac  IC  519.07-lB). 
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that  altitude.  The  theoretical  value  at  180  km  is  3 to  4 times  larger  than 
the  radar  data.  If  our  input  particle  spectrum  (Figure  2-6)  is  a reasonable 
representation  of  the  precipitation  into  che  arc  system,  then  we  are  unable 
to  account  for  the  difference.* 

3914-A  EMISSION 


In  Figure  2-10  the  zenith  intensity  at  3914A  measured  by  the  rocket- 
borne  photometer  (Reference  10)  is  compared  with  the  ARCTIC  code  calculations 
and  with  the  intensity  inferred  from  the  SRI  radar  data.  The  latter  was 
obtained  by  first  computing  a steady  state  volume  emission  rate  and  then 
integrating  it  over  altitude.  The  volume  emission  rate,  is  propor- 

tional to  the  energy  deposition  rate  which,  in  steady  state,  is  proportional 
2 

to  n . Therefore, 


'3914  ARTIC 


^"e  ^ARCTIC 


where  the  subscripts  "ARCTIC"  and  "SRI"  refer  to  values  obtained  from  the 
ARCTIC  code  and  from  the  SRI  radar,  respectively.  An  altitude  plot  of 

o 

In  comparing  the  3914-A  rocket  data  (ascent)  with  the  code  results, 
it  should  be  borne  in  mind  that  the  rocket  did  not  enter  the  main  arc  until 
it  was  at  about  125-km  altitude  (see  Figure  2-3).  When  it  was  below  this 
altitude,  only  part  of  the  photometer  sight  path  intersected  the  main  arc 

* An  overestimate  of  the  input  flux  for  electron  energies  between  0.1  and 
0.7  keV  cannot  account  for  it.  The  input  flux  in  this  energy  range  is 
2.4x10^^  eV  cm"2  sec"^.  The  energy  deposition  rate  required  to  produce, 
in  steady  state,  the  SRI  and  ARCTIC  code  profiles  of  ng  above  120-km 
altitude  are,  respectively,  6.8  xio^^  and  2.0  xio^^  ev  cm"^  sec"l.  The 
difference  is  1.32 x 10^^  eV  cm"^  sec'^  and  this  is  55  times  larger  than 
the  assumed  incident  energy  flux  in  the  0.1  to  0.7  keV  range. 
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Figure  2-10.  Zenith  radiance  at  3914A  (Sergeant  Hydac  IC519.07-1B) ; I 

data  and  theory  compared.  I 


complex.  The  calculations,  however,  are  based  on  a vertical  sight  path 
entirely  within  the  arc.  Therefore,  until  the  rocket  entered  the  arc,  one 
might  expect  the  measured  intensity  to  be  somewhat  lower  than  the  calculated 
values.  The  comparison  shows,  however,  that  the  measured  intensity  is 
almost  uniformly  lower  (factor  3 to  4)  than  that  calculated  over  the  height 
interval  from  90  to  175  km.  This  discrepancy  is  consistent  with  earlier 
comparisons,  especially  those  for  the  '73  and  '74  Multi  events  reported 
in  References  2 and  4.  The  reason  for  this  persistent  difference  between 
the  photometer  data  and  our  calculations,  based  on  the  particle  data,  is 
not  known . 

O 

The  3914A  intensity  profile  inferred  from  the  radar  data,  shown 
in  Figure  2-10,  is  also  larger  than  the  rocket  data  below  about  120  km. 
However,  above  about  125-km  altitude  it  is  smaller,  thus  tending  to  sub- 
stantiate the  notion  that  the  energy  deposition  rate  inferred  from  the 
radar  data  shown  in  Figure  2-9  may  be  too  small  in  this  altitude  region. 

5.3-ym  EMISSION 


Figure  2-11  shows  a comparison  between  the  rocket  data  (Reference 
10)  at  5.3  ym  and  the  code  calculation  of  NO  chemiluminescence.  Both  the 
radiometer  and  the  CVF  spectrometer  data  are  shown.  The  latter  was 
obtained  by  integration  over  the  spectral  features  from  5.0  to  5.6  ym. 

The  integrated  CVF  data  are  seen  to  be  larger  than  the  radiometer  data  by  a 
factor  of  about  2 on  rocket  descent  and  somewhat  less  than  that  on  ascent. 
Part  of  this  difference  is  probably  attributable  to  the  fact  that  the 


integrated  CVF  data  covers  a band  of  width  0.6  ym,  whereas  the  full  width 
at  half  maximum  of  the  radiometer  response  (Reference  11)  is  only  about 
0.3  ym.  The  solid  curve  in  Figure  2-11  is  the  calculated  contribution  to 
the  emission  from  NO  chemiluminescence.  It  is  the  total  calculated  band 
emission  reduced  by  a factor  of  0.23  (as  explained  in  Reference  2,  p.  65)  to 
to  allow  for  the  radiometer  relative  response  function. 

28 


Ascent - 


-Descent 


» **  X 


®N  • \ 

o\  \ 


\ \ 

\ 

\ s\ 


ARCTIC  (Chemi  - 100 ^sec) 

ARCTIC  Chemi  + Descent  Radiometer  Data 

ARCTIC  Chemi  + Descent  CVF  Data 

o o o Radiometer  Data  | 

« « « CVF  Data  (Integrated) 


80  100  120  140  160  180  160  140  120  100  80 

Rocket  Altitude  (km) 

Figure  2-11.  Calculated  and  observed  zenith  radiance  in  5.3-vim  band 
(Sergeant  Hydac  IC  519.07-lB). 


■1 


The  calculated  chemiliuninescence  (solid  curve  in  Figure  2-11)  is 
seen  to  be  smaller  than  the  ascent  data.  However,  if  we  add  to  it  the  radiance 
measured  by  the  radiometer  on  rocket  descent,  the  resulting  curve  (dashed) 
is  in  reasonable  agreement  with  the  ascent  radiometer  data.  Similarly,  if 
i we  add  to  the  chemiluminescence  the  radiance  measured  by  the  CVF  spectrometer 

on  descent,  the  resulting  curve  (dot  dash)  is  in  reasonable  agreement  with 
the  ascent  CVF  data,  especially  below  about  130  km.*  The  rationale  for 
doing  this  is  the  assumption  that  the  descent  data,  taken  in  a region  of 
low  particle  precipitation,  is  attributable  to  NO  emission  excited  by  the 
atom  interchange  reaction  N0(v=0)  + 0' -►  N0'(v=l)  +0  and  by  earthshine 
scatter  from  NO  that  operate  equally  on  ascent  through  the  arc  and  on 
descent  outside  the  arc.  Thus,  we  assume  that  the  proper  amount  to  add  to 
the  calculated  chemiluminescence  on  ascent  to  allow  for  these  additional 
excitation  mechanisms  is  just  the  measured  descent  radiance.  The  NO 
profile  required  to  match  the  descent  radiometer  data  is  shown  in  Figure  5-1. 


t 


I 

I 
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It  should  be  emphasized  that  the  chemiluminescent  contribution  to 
the  NO  emission  is  based  on  certain  assumptions  that  have  not  yet  been 
verified  experimentally.  The  principal  assumption  is  that,  following 
chemical  formation  of  NO,  the  vibrational  states  are  populated  with  equal 
probability  up  to  the  maximum  allowed  by  energy  conservation.^  This  leads 
to  the  emission  spectrum  shown  in  Figure  2-9  of  Reference  4.  A different 
assumption  would  have  led  to  a different  spectrum  and  to  a different  fraction 
of  the  energy  emitted  in  the  radiometer  pass  band  and  in  the  5.0  to  5.6-pm 
wavelength  interval. 


* For  comparison  with  the  CVF  data,  we  have  taken  the  chemiluminescent  con- 
tribution equal  to  0.44  of  the  total  calculated  band  emission.  This  is 
because  44  percent  of  the  energy  in  the  calculated  chemiluminescent 
spectnam  (Figure  2-9  of  Reference  4)  is  contained  in  the  wavelength  band 
from  5.0  to  5.6  ym. 

^ In  our  model,  most  of  the  NO  chemiluminescence  arises  from  the  reaction 
N(^D)  +02-*N0+  + 0, 


tkitt.. 
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The  CVF  spectra  provide  some  basis  for  understanding  the  emission 
mechanisms  and  radiating  species  involved,  although  their  usefulness  in  this 
regard  is  limited  by  the  instrument  resolution  and  the  long  wavelength  cut- 
off at  about  5.64  pm.  Assuming  that  the  (1,0)  band  of  NO  is  a principal 
contributor,  we  have  calculated  what  the  corresponding  CVF  spectrum  should 
look  like  for  a rotational  temperature  of  250  “K.  The  results  are  shown 
in  Figure  2-12.  Three  different  spectra  are  shown,  each  representing  a 
different  instrument  resolution.  For  the  instrument  used  on  rocket 
IC  519.07-lB,  the  resolution  near  5.3  ym  is  about  0.13  ym  (Reference  12). 
However,  our  theoretical  curves  are  based  on  a flat  response  over  the 
resolution  element  whereas,  in  practice,  it  is  more  nearly  triangular,  or 
at  least  gaussian.  Therefore,  with  respect  to  our  calculated  spectra,  the 
lower  one  in  Figure  2-12,  corresponding  to  0.12-ym  resolution,  is  probably  the 
most  applicable. 

Figure  2-13  compares  this  theoretical  spectrum  with  a typical  CVF 
feature.  The  two  spectra  have  been  normalized  so  that  the  magnitude  of 
their  peaks  coincide.  Although  the  observed  feature  is  broader  than  the 
theoretical  one,  the  positions  of  the  CVF  peaks  are  in  striking  agreement 
with  the  R-  and  P-branch  peaks  of  the  (1,0)  band  of  NO. 

One  might  be  tempted  to  attribute  the  broader  CVF  features  to 
contributions  from  higher  vibrational  transitions  of  the  Av=l  sequence. 
However,  although  this  could  account  for  broadening  on  the  long  wavelength 
side  of  the  feature,  it  could  not  account  for  broadening  on  the  short 
wavelength  side  unless  the  rotational  temperature  is  effectively  much 
higher  than  what  we  have  assumed.  In  fact,  it  is  likely  that  the  emission 
at  5.0  ym,  which  appears  persistently  in  CVF  data  from  all  ICECAP  rockets, 
arises  from  some  source  other  than  NO.* 


* The  suggestion  has  been  made  that  CO2  may  be  a principal  contributor 
near  5.0  ym.  However,  the  altitude  dependence  of  the  emission  does  not 
support  this  hypothesis. 
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Figure  2-13.  Typical  CVF  spectrometer  feature  near  5.3  pm  (from  rocket 
IC  519.07-lB)  compared  with  the  theoretical  (1,0)  band 
spectrum  of  NO  at  250  'K. 
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The  observed  CVF  features  are  not  consistent  with  the  calculated 
NO  chemiluminescent  spectrum  (based  on  v = 0 to  18  equally  populated)  shown 
in  Reference  4 (Figure  2-9b) . This  spectrum  exhibits  considerable  emission 
at  wavelengths  longer  than  the  spectrometer  cut-off  at  5.64  yra.  But  the 
observed  features  usually  drop  below  noise,  even  at  5.60  ym,  indicating 
little  or  no  signal  at  5.64  ym. 

By  contrast,  as  we  have  seen  above,  the  data  are  much  more  con- 
sistent with  the  assumption  that  only  the  (1,0)  band  of  NO  is  principally 
involved.  This  supports  the  atom  interchange + earthshine  scatter  hypothesis 
for  which  only  N0(v=  1)  would  be  appreciably  excited. 

4.3-ym  EMISSION 

Figure  2-14  shows  the  measured  zenith  radiance  near  4.3  ym.  It 
includes  the  radiometer  data  and  also  the  CVF  spectrometer  data  integrated 
in  wavelength  over  the 4.3-ym  features.  Again,  the  wavelength  intervals 
covered  in  Figure  2-14  by  the  two  instruments  are  not  identical.  The  half 
power  points  on  the  radiometer  response  (Reference  11)  lie  between  4.21  and 
4,31  ym,  whereas  the  CVF  features  were  integrated  from  4.08  to  4.32  ym. 

Also  shown  in  Figure  2-14  are  several  theoretical  results.  The 
dashed  curves  give  the  calculated  zenith  radiance  from  CO2  vibraluminescence 
after  100  and  1000  sec  of  arc  bombardment.  The  dotted  curve  is  Kumer's 
estimate  (Reference  13)  of  the  ambient  nighttime  CO2  radiance  from  thermal 
collisions  and  earthshine  scatter.  The  dot-dash  curve  is  an  estimate  of  the 
CO2  vibraluminescence  due  to  drizzle  bombardment  of  the  ascent  rocket  air- 
space during  the  half  hour  prior  to  rocket  launch.  It  is  based  on  an  ARCTIC 

code  calculation  of  the  drizzle  environment  assumed  for  the  '74  Multi  event 

_2 

(Reference  4)  and  results  from  an  energy  input  flux  of  about  4 ergs  cm 
sec"^.  This  is  reasonably  consistent  with  the  measured  flux  shown  in 


Figure  2-4,  although  it  may  be  a little  on  the  low  side.  The  solid  curve  is 
the  sum  of  the  contributions  from  the  arc  (100  sec),  the  drizzle,  and  the 
ambient  CO2  radiance.  The  ARCTIC  code  results  (dashed  and  dot-dash  curves) 
have  been  corrected  for  the  radiometer  relative  response  function  so  as  to 
make  them  comparable  with  the  radiometer  data. 

The  solid  curve  is  seen  to  be  in  good  agreement  with  the  radiometer 
data  between  about  90  and  105  km.  However,  at  higher  altitudes  the  calcu- 
lated vibraluminescence  falls  off  faster  than  does  the  ascent  radiometer 
data  so  that  by  110  km  they  differ  by  a factor  of  2. 

It  should  be  realized  that  detailed  agreement  between  auroral 
data  at  4.3  pm  and  a calculation  of  CO^  vibraluminescence,  such  as  we  have 
done,  is  not  generally  to  be  expected.  The  reason  for  this  is  that  energy 
deposited  in  vibrational  excitation  of  at  any  given  time  is  (partially) 
transferred  to  CO^  and  radiated  at  4.3  pm  with  an  e-folding  time  ranging 
from  about  8 minutes  at  100  km  to  over  one  half  hour  at  120  km.  Therefore, 
the  CO2  emission  at  any  given  time  is  the  integrated  contribution,  due  to 
vibrational  excitation  of  N2,  over  past  times.  An  accurate  calculation 
would  require  a knowledge  of  the  spectral  flux  of  precipitating  electrons 
for  minutes  to  hours  (depending  on  altitude)  prior  to  the  measurement  time. 

In  addition,  as  pointed  out  by  Kumer  (Reference  13),  transport  effects  due 
to  winds  and  diffusion  can  appreciably  modify  the  local  concentration  of 
N2(v  = 1),  and  hence  the  CO2  emission,  in  the  time  scales  involved.  These 
effects  should  also  be  included.  Detailed  information  of  this  nature  has 
yet  to  be  available  in  any  auroral  experiment. 

In  the  calculations  presented  in  Figure  2-14,  we  used  the  measured 
arc  flux  and  guessed  that  it  bombarded  the  rocket  airspace  for  100  sec. 

This  seems  reasonable  on  the  basis  of  the  arc  morphology  inferred  from  the 
Lockheed  photometer  data,  but  the  magnitude  of  the  emission  is  sensitive 
to  the  actual  bombardment  time  which  could  have  been  at  least  a factor  of  2, 


36 


more  or  less,  than  what  we  assumed.  The  predosing  from  the  drizzle  is  also 
a guess,  albeit  aided  by  the  ground  photometer  data.  Furthermore,  no  cor- 
rection for  transport  effects  has  been  made.  Considering  all  these  un- 
certainties and  assumptions,  the  agreement  with  the  data  should  be  considered 
good.  We  feel  that  the  CO^  vibraluminescence  can,  in  the  main,  account 
for  the  ascent  data  at  4.3  ym.  It  is  likely,  however,  that  on  ascent 
through  the  arc  complex,  the  instruments  received  additional  contributions 
at  4.3  ym  from  NO^  chemiluminescence.  This  possibility  is  discussed  below. 


In  Figure  2-14  the  descent  data  near  4.3  ym  are  also  shown.  It 
should  be  noted  that  although  the  rocket  descended  in  a region  of  low 
particle  precipitation,  the  radiance  measured  near  4.3  ym  is  almost  as 
large  as  that  on  rocket  ascent  through  the  arc.  This  is  further  evidence 
in  support  of  CO2  vibraluminescence  as  the  main  mechanism  operating  near 
4.3  ym.  For  if  the  radiance  were  mainly  from  NO^  chemiluminescence,  then 
the  data  on  descent  through  the  drizzle  should  have  been  considerably  less 
than  that  on  ascent  through  the  arc  because  the  time  constant  for  decay 
of  NO^  chemiluminescence  is  short  (a  few  seconds) . On  the  other  hand, 
the  rocket  descent  airspace  had  been  bombarded  by  the  arc  flux  about  12 
minutes  earlier  as  the  arc  system  moved  south  through  it.  The  decay  of  the 
CO2  vibraluminescence  is  sufficiently  slow  so  that  the  emission  seen  on 
descent  could  result  from  this  earlier  bombardment. 


Figure  2-15  shows  the  calculated  zenith  intensity  at  4.3  ym  from 
NO^  chemiluminescence  in  the  arc.  The  upper  dashed  curve  is  an  integral 
over  altitude  of  the  volume  emission  rate  shown  in  Figure  2-8.  It  is  for  the 
entire  band,  uncorrected  for  instrument  response,  and  it  does  not  include 


collisional  quenching.  Clearly,  this  uncorrected  intensity  is  much  larger 
than  that  observed  by  the  radiometer.  The  middle  dashed  curve  in  Figure  2-15 
incorporates  a correction  factor  of  0.16  to  allow  for  the  radiometer  relative 


response  function  (see  Reference  2, 


p.  72). 


This  correction  factor  is  based 
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on  the  assumption  that  the  chemiluminescent  spectrum  shown  in  Figure  2-16 
(Reference  4)  is  applicable.  No  laboratory  data  exist  to  support  it. 

The  lowest  curve  in  Figure  2-15  (dot-dash)  incorporates  an  addi- 
tional correction  for  collisional  quenching.  To  our  knowledge,  quenching 
rates  for  the  vibrational  states  of  NO^  have  not  been  measured.  However, 
if  we  assume  that  vibrational  quenching  of  NO*  proceeds  mainly  by  v-v 
transfer  to  N2  with  rate  constant  k,  then  we  can  allow,  approximately, 
for  it  by  reducing  the  volume  emission  rate  by  the  factor  Aj^Q/(Aj^Q+k[N2]) 
where  A^^  is  the  Einstein  transition  rate  for  the  fundamental  band  of 

NO*.  The  dot-dash  curve  is  based  on  values  for  A and  k of  13.4  sec’^ 

-9  3 -1 

and  1 X 10  cm  sec  , respectively.  When  the  dot-dash  curve  is  added 
to  the  CO^  emission  (solid  curve  in  Figure  2-14),  the  solid  curve  in  Figure 
2-15  results.  The  radiometer  data  are  shown  for  comparison- 

Unfortunately,  as  explained  above,  there  are  too  many  uncertainties 

inherent  in  the  CO2  calculation,  as  well  as  in  the  photon  efficiency 

spectrum  of  NO  chemiluminescence,  to  say  if  the  apparent  agreement  with 

the  4.3-ym  data  shown  in  Figure  2-15  has  any  significance.  The  only  thing 

we  can  say  for  sure  is  that  if  the  assumptions  inherent  in  our  CO2  and 

NO*  calculations  are  reasonably  good,  then  a (lumped)  quenching  rate 

+ -93-1 

constant  for  the  vibrational  states  of  NO  as  large  as  1 x 10  cm  sec 
is  needed.  A value  of  this  magnitude  is  not  unusual  for  an  ion-molecule 
reaction. 

2,7-ym  EMISSION 


No  data  at  2.7  ym  are  available  from  rocket  IC  519.07-lB  with 
which  to  compare  our  calculated  NO  chemiluminescence.  The  CVF  spectrometer 
data  in  this  spectral  region  appear  to  be  primarily  noise.  This  is  consistent 
with  our  calculated  intensity  at  2.7  ym  that  is  below  the  instrument  noise 
level  of  2.3  MR(ym)  To  this  extent,  therefore,  our  calculated  NO  chemi- 
luminescence is  not  in  disagreement  with  the  CVF  spectrometer  observations. 
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SUMMARY  AND  CONCLUSIONS 


Results  and  conclusions  from  our  study  of  the  12  March  1975 
auroral  event  can  be  summarized  as  follows; 

O 

(1)  The  calculated  3914-A  intensity  profile  for  the  arc  system, 
based  on  the  measured  particle  flux,  is  a factor  of  3 to  4 times  larger 
than  that  measured  by  the  rocket-borne  photometer.  This  result  is  con- 
sistent with  corresponding  comparisons  made  for  other  auroral  events, 
although  the  reason  for  it  is  not  known. 

(2)  The  emission  intensity  observed  near  5.3  pm  on  rocket  descent 
through  a low  level  drizzle  precipitation  region  can  be  accounted  for  on 
the  basis  of  atom  interchange  + earthshine  scatter  excitation  of  NO  for 
a suitably  chosen  enhanced  NO  profile.  A combination  of  this  descent 
emission  and  the  calculated  NO  chemiluminescence  in  the  arc  is  in  good 
agreement  with  the  data  measured  on  rocket  ascent  through  the  arc  complex. 
The  peaks  in  the  CVF  spectra  near  5.3  pm  correspond  closely  to  the  R and 

P branches  of  the  (1,0)  fundamental  band  of  NO.  Thus,  the  data  are 
generally  consistent  with  NO  emission,  excited  by  a combination  of  atom 
interchange  and  earthshine  scatter,  as  the  principal  source  of  radiation 
near  5.3  pm. 

(3)  The  ascent  and  descent  rocket  data  near  4.3  pm  are  consistent 
with  CO^  vibraluminescence  as  the  dominant  emission  mechanism.  However, 
detailed  agreement  with  theory  cannot  be  expected  unless  accurate  account 
is  taken  of  predosing  and  the  transport  of  vibrationally  excited  N^  by 
winds  and  diffusion. 

(4)  In  order  to  suppress  NO*  chemiluminescence,  that  would  other- 
wise dominate  the  emission  near  4.3  pm,  a quenching  rate  coefficient  (by 

-9  3 -1 

N^)  as  large  as  1 x 10  cm  sec  is  needed. 

(5)  The  CVF  spectrometer  data  near  2.7  pm,  consisting  of  noise,  is 
not  inconsistent  with  our  calculated  emission  from  NO  chemiluminescence 
that  also  provides  for  signal  levels  below  the  instrument  noise. 


SECTION  3 

AURORAL  EVENT  OF  10  MARCH  1975 
(Nike  Hydac  IC  507.11-2A) 


PRELIMINARY 

On  10  March  1975,  at  0912:20  UT  (2312:20  local  time),  rocket 
IC  507.11-2A  was  launched  from  Poker  Flat  into  a very  active  auroral  breakup, 
reaching  an  apogee  of  152.3  km  at  200  sec  after  launch.  The  region  traversed 
by  the  rocket  was  the  central  intense  part  of  a westward  traveling  surge 
(Reference  5) . The  scanning  photometer  data  of  Sears  (Reference  9)  showed 
that  the  region  was  very  broad,  extending  30  degrees  on  either  side  of  the 
zenith  at  Chatanika. 

The  rocket-borne  instrumentation  included  three  vertical -viewing, 
helium-cooled  radiometers  operating  in  the  5.6-,  9.6-,  and  7-12-  pm  regions. 

No  particle  detectors  or  other  optical  instrumentation  were  on  board. 

In  addition  to  the  Lockheed  scanning  photometer  data  (Reference  9) , 
other  ground  support  data  for  this  event  are  available  from  the  SRI  incoherent 
backscatter  radar  (Reference  6) . Camera  malfunction  precluded  any  scanning 
photometer  data  from  Ft.  Yukon  or  Ester  Dome. 

Since  no  measured  particle  flux  data  are  available  for  this  event, 
and  only  very  limited  optical/IR  data,  no  attempt  was  made  to  make  a detailed 
ARCTIC  code  calculation.  Rather,  we  have  attempted  to  scale  from  previous 
calculations  to  infer  what  the  zenith  intensity  due  to  NO  chemiluminescence 
should  have  been,  and  to  compare  the  results  with  the  radiometer  data. 
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The  total  energy  deposition  rate  from  precipitating  electrons  has 

been  inferred  from  the  Lockheed  photometer  data  (Reference  9)  and  the  SRI 

radar  data  (Reference  6) . Both  sets  of  data  show  the  deposition  rate  peaking 

at  about  0915  UT,  just  after  rocket  launch.  However,  the  photometer  data 

-2  -1 

shows  a peak  of  over  200  ergs  cm  sec  , while  the  radar  data  indicate 
-2  -1 

about  62  ergs  cm  sec  , over  a factor  of  3 less. 

Since  the  radar  data  provide  an  altitude  profile  of  the  electron 
density,  n^,  we  have  used  it  to  infer  an  energy  deposition  rate  profile,  e. 

In  particular,  in  steady  state, 

e = nan^^  (3-1) 

where  n is  eV/ion  pair,  and  a is  the  effective  electron  recombination 
coefficient.  Figure  3-1  shows  the  electron  density  profile  (Reference  6) 
measured  for  the  time  period  0914:22-0920:58.  Since  it  is  fairly  similar 
to  the  one  calculated  for  the  *73  Multi  arc,  we  have  used  our  results  from 
that  calculation  to  provide  profiles  for  n and  a,  the  latter  of  which 
is  shown  in  Figure  3-1.  The  value  of  n varies  slowly  from  about  35.1  at 
90  km  to  33.4  at  175  km.  The  resulting  profile  of  e is  also  shown  in 
Figure  3-1 , 

COMPARISONS  WITH  DATA  IN  5.3-Mm  BAND 

To  compute  the  NO  chemiluminescence  from  the  deposition  rate 
profile  shown  in  Figure  3-1,  we  have  adopted  the  5.3-vim  photon  efficiency 
(photons  per  eV  deposited)  calculated  for  the  *73  Multi  event  (rocket 
A18. 205-1)  for  a bombardment  time  of  80  sec.  This  is  shown  in  Figure  3-2, 
along  with  the  resulting  volume  emission  rate  profile. 
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Figure  3-2.  Calculated  volume  emission  rate  near  5.3  ym  from  NO 
chemiluminescence  for  auroral  event  of  10  March  1975. 
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Figure  3-3  shows  a comparison  between  the  intensity  measured  by 
the  radiometer  on  rocket  ascent  and  our  computed  intensity  from  NO  chemi- 
luminescence. (The  descent  data  are  similar,  at  least  above  110  km) . The 
effective  radiometer  passband  covers  the  wavelength  interval  from  5.27  to 
5.88  ym  (Reference  14).  The  theoretical  curve  in  Figure  3-3  is  the  altitude- 
integrated  volume  emission  rate  (Figure  3-2)  with  a correction  factor  of  0.65 
applied  to  allow  for  the  fraction  of  the  chemiluminescent  spectrum  (Reference 
2,  Figure  A-1)  that  falls  in  the  radiometer  passband. 

The  agreement  between  the  two  curves  in  Figure  3-3  is  obviously 
poor,  especially  above  100  km  where  the  theoretical  curve  falls  off  with 
increasing  altitude  much  faster  than  does  the  data.  One  might  argue  that 
this  may  be  due  to  an  underestimate  of  the  energy  deposition  at  higher 
altitude^,  based  on  the  radar  data,  as  appeared  to  be  the  case  for  the 
12  March  1975  event  (see  Section  2,  especially  Figure  2-9).  However,  our 
revised  ARCTIC  code  calculation  for  the  '73  Multi  event  gives,  for  80  sec 
of  arc  bombardment,  an  NO  emission  profile  similar  in  shape  to  that  shown 
in  Figure  3-3  and  only  a factor  of  2 larger.  It  is  unlikely,  therefore, 
that  had  a measured  particle  spectrum  been  available  to  us,  and  a full  code 
calculation  been  performed,  results  much  different  from  those  shown  in 
Figure  3-3  would  have  been  obtained.  In  fact,  all  previous  runs  with  the 
ARCTIC  code  have  shown  that,  although  the  magnitude  of  the  chemiluminescent 
intensity  changes  as  the  input  spectral  flux  is  changed,  the  shape  of  the 
altitude  profile  remains  approximately  the  same  as  shown  in  Figure  3-3. 

It  does  not  seem  possible  to  match  the  shape  of  the  radiometer  curves  in 
Figure  3-3  with  NO  chemiluminescence  alone. 


It  is  tempting,  as  in  previous  comparisons  with  data  near  5.3  ym,  j. 


to  attribute  the  difference  to  excitation  of  NO  by  a combination  of  the 
atom  interchange  reaction,  N0(v=0)  + 0' -►N0''(v=l)  +0,  and  earthshine  scatter, 
operating  on  an  aurorally  enhanced  NO  profile.  However,  if  this  is  indeed 
the  explanation,  then  NO  concentrations  4 to  5 times  larger  than  the 
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already  enhanced  values  that  seem  to  have  been  present  under  certain  other  | 

auroral  conditions  would  be  required.  This  follows  from  a comparison  with  | 

the  radiometer  data  obtained  on  descent  of  the  '73  and  '75  Multi  rockets  I 

through  drizzle  regions  (shown  by  the  dashed  curves  in  Figure  3-3) . These  I 

data  have  previously  been  ascribed  to  NO  emission  excited  by  earthshine  | 

scatter  and  atom  interchange,  and  are  seen  to  be  4 to  5 times  smaller  than  1 

the  present  radiometer  data.*  In  fact,  if  we  ascribe  all  of  the  present  j 

radiometer  data  to  these  mechanisms,  the  inferred  NO  profile  is  shown  in  I 

Figure  5-1  of  Section  5.  It  is,  indeed,  seen  to  be  considerably  enhanced  I 

over  the  other  enhanced  profiles.  | 

] 

] 

Although  intense  precipitation  occurred  during  the  rocket  flight,  I 

its  duration  in  the  rocket  airspace  prior  to  launch  was  too  short  for  any  1 

appreciable  buildup  of  the  NO.  Probably  of  more  significance,  however,  i 

is  the  predosing  received  over  longer  time  periods  prior  to  rocket  launch. 

We  have  estimated  the  relative  predosing  for  the  present  event  (of  10  March  j 

1975)  and  for  the  '75  Multi  event  (2  days  later)  using  the  time  histories 

of  energy  deposit  by  particle  precipitation  given  by  Sears  (Reference  9)  • 

and  Vondrak  (Reference  6)  . From  the  optical  data  (Reference  9) , the  ratio 

of  the  integrated  energy  deposition  over  a time  span  of  67  min  prior  to 

rocket  launch  on  10  March  1975,  to  that  on  12  March  1975,  is  about  3.2. 

From  the  radar  data  (Reference  6),  the  corresponding  ratio,  over  a time 
span  of  138  min,  is  4.3. 

These  numbers  suggest  that  the  integrated  dosage  prior  to  the 
10  March  event  may  have  been  sufficiently  greater  than  that  prior  to  the 
12  March  event  to  produce  the  additional  NO  concentration  required. 

Longer  integration  times  are  needed,  however,  to  say  with  more  assurance  1 

whether  or  not  this  is  likely  to  be  the  case.  ; 



* Although  the  radiometer  used  for  the  10  March  1975  event  was  not  identical 
to  those  used  in  the  other  two  experiments,  we  feel  that  the  difference 
is  not  enough  to  be  very  significant. 
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SUMMARY  AND  COMMENTS 


A 


Based  on  the  electron  density  profile  measured  by  the  SRI  back- 
scatter  radar,  we  have  inferred  an  energy  deposition  rate  profile  for  the 
auroral  event  of  10  March  1975  and,  together  with  results  from  previous 
ARCTIC  code  calculations,  have  estimated  the  intensity  profile  of  NO 
chemiluminescence  near  5.3  ym.  The  results  show  that  chemiluminescence 
cannot  account  for  the  radiometer  data,  especially  at  altitudes  above  100  km. 


It  is  possible  to  account  for  the  results  on  the  basis  of  NO 
emission  excited  by  earthshine  scatter  and  the  atom  interchange  reaction, 
but  only  by  invoking  an  NO  profile  that  is  considerably  enhanced  over 
already  enhanced  profiles  inferred  from  other  auroral  events.  A limited 
survey  of  data  on  predosing  suggests  that  the  integrated  energy  deposited 
in  the  region  prior  to  rocket  launch  may,  at  least  qualitatively,  account 
for  the  presence  of  larger  NO  values.  However,  for  this  very  intense 
event,  one  would  have  expected  NO  chemiluminescence  near  5.3  ym  to 
dominate  other  mechanisms.  It  is  rather  disconcerting  that  this  does  not 
appear  to  have  been  the  case. 
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SECTION  4 

CODE  CALCULATION  OF  THE  ELECTRON  TEMPERATURE 
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PRELIMINARY 


In  comparing  ARCTIC  code  calculations  of  electron  temperature, 
performed  prior  to  the  present  effort,  with  observational  data  for  the 
auroral  environment,  it  became  clear  that  the  code  results  were  too  large. 
In  fact,  at  low  altitudes  (~100  km)  in  strong  arcs  (such  as  the  ’74  Multi 
event)  our  temperature  values  were  high  by  as  much  as  a factor  of  5.  Such 
errors  can  impact  severely  upon  the  accuracy  of  calculated  species  concen- 
trations, especially  those  of  the  ions  and  electrons.  Considerable  effort 
was  devoted  to  tracking  down  the  source  of  the  error  and  in  providing  a 

t 

computational  method  that  we  now  believe  is  accurate.  It  incorporates 
several  features  that  have  not  hitherto  been  reported  in  the  literature. 
This  section,  along  with  Appendices  A and  B,  is  devoted  to  a description 
of  this  work. 

/ 

HEATING  AND  COOLING  RATES  / 

The  basic  equation  used  in  the  ARCTIC  code  for  determining  the 
electron  temperature  is  given  by  (Reference  1) 

ar  ■ 

where  H and  C are,  respectively,  the  heating  and  cooling  rates  per 

electron.  The  ARCTIC  code  includes  the  following  basic  cooling  terms: 

1)  excitation  of  0*(^D),  0'"(^P),  2)  excitation  of  0(^0) , 0(^S),  N(^D), 

2 

N{  P) , 3)  excitation  of  vibrational  states  of  N2,  4)  excitation  of 


.ji 
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vibrational  states  of  O2,  5)  excitation  of  fine  structure  levels  of  0(^P) , 
6)  rotational  excitation  of  O2.  The  details  of  these  cooling  terras  can 

be  found  in  Reference  1. 

In  addition  to  electron  cooling  by  these  processes,  we  have 

recently  added  cooling  from  elastic  scattering  by  N , O2,  0 and  positive 

-3  ^ -1 

ions.  Explicit  expressions  for  the  rates  (eV  era  sec  ) are  given,  from 
Reference  16,  by: 


= 1.77X10'^^  (2/3k)  [N,]  (1-1.21  X 10"'^  T )T  (T  -T  ) 

Z.  6 6 6 ri 

(4-2) 

= 1.21x10'^®  (2/3k)  [0-]  (1+3.6  X lO'^T  *^)T  ^(T  -T  ) 

2-'^  e e en 

(4-3) 

= 7.9x10"^^  (2/3k)  [0]Tg*^(Tg-Tj^)  (1+5.7  X lo'^^Tg) 

(4-4) 

c. 

10ns 

7 T - T 

= 5x10'^  (2/3k)  -^-37!  ([0  ] + .5[02'']  + .53[N0^]) 

^e 

(4-5) 

where  2/3k  = 7.7x10^  '’K/eV. 
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where  d«^/dt  is  the  rate  energy  is  deposited  in  the  electron  gas.  It  is 
determined  from  the  expression 


J <ti(E)L^(E.Tg)dE 


(4-8) 


o 

where  <{»(E)  represents  the  net  steady  state  flux  (both  primaries  and 

secondaries)  and  is  the  electron  loss  function  (see  Appendix  A).  In 

addition,  heating  of  the  electron  gas  also  occurs  due  to  deexcitation  of 

□■^(^D),  O'^(^P),  0(^D),  0(^S),  N(^D),  and  N(^P)  . However,  except  for 
2 

deexcitation  of  N(  D)  at  high  altitudes  (~200  km),  these  heating  terms 
involve  only  small  energy  transfer.  Furthermore,  not  all  of  the  energy  of 
the  ejected  electron  ends  up  in  the  ambient  electron  gas.  Appendix  B gives 
a description  of  how  we  partition  the  ejected  electron  energy. 


OLD  SOLUTION  TO  EQUATION  4-1 


In  order  to  solve  Equation  4-1  for  T^,  we  have  to  solve  Equations 
4-7  and  4-8  for  the  heating  terms.  If  this  is  done  every  numerical  time  step 
in  the  chemistry  code,  it  becomes  prohibitively  expensive.  Therefore, 
approximate  solutions  are  usually  adopted  for  Equations  4-7  and  4-8.  In 
Appendix  A we  show  that  the  usual  approximation  made  (see  also  Reference  17) 
is 

dS 

« d^j/dt  (4-9) 

where  f^  is  the  fractional  ionization,  6 is  a constant,  and  d<^^/dt  is 
the  total  energy  deposition  rate.  Since  the  fractional  ionization  is  given 
by  n^/A,  where  A is  the  total  particle  concentration.  Equation  4-7 
becomes 

H = 7.7x10^  :i-d<f^/dt  . (4-10) 
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Thus,  with  approximation  4-9,  the  heating  rate  per  electron  is  independent 
of  the  electron  density. 

This  is  the  assumption  previously  made  in  the  ARCTIC  code  and, 
as  far  as  we  are  aware,  by  all  other  workers  in  the  field.  At  a given 
altitude,  the  code  would  perform  a detailed  calculation  of  6 for  a 
given  electron  density  and  would  assume  the  validity  of  Equation  4-10  during 
the  temperature  and  chemistry  calculation  while  the  electron  density  varied. 
Since  the  cooling  rate  is  also  independent  of  the  electron  density,  the 
electron  temperature  was  then  also  independent  of  the  electron  density. 

I'  -Jer  this  approximation,  the  solution  to  Equation  4-1  can  be  represented 
schematically  as  shown  in  Figure  4-1.  We  see  that  it  exhibits  a monotonic 
rise  following  onset  of  particle  bombardment  until  it  reaches  a constant 
steady  state  value  can  be  found  by  setting  Equation  4-1  to 


-Steady  State  Temperature  T 


Electron  Gas  Temperature  T^ 
y — Neutral  Gas  Temperature 


Figure  4-1.  Behavior  of  Tg  and  Tgj  under  approximations  4-9  and  4-10. 


zero.  Often  the  monotonic  rise  can  be  represented  by  an  exponential  rise 
with  a time  constant  given  by 


T =T 
' e ss. 


(4-11) 


In  the  ARCTIC  code,  detailed  account  was  taken  of  the  rise  portion 
of  Figure  4-1.  However,  once  the  electron  temperature  reached  we  set 

it  equal  to  for  all  subsequent  times.  As  mentioned  above,  this  pro- 

cedure has  resulted  in  gross  errors. 


NEW  SOLUTION  TO  EQUATION  4-1 


A detailed  study  of  the  validity  of  Equations  4-9  and  4-10  has 

been  made.  This  work  is  reported  in  Appendix  A.  Basically,  the  results  of 

-9 

this  study  show  that  Equations  4-9  and  4-10  hold  for  f < 10  . For 

-9  -6  ® 

10  i f £ 10  , the  partition  of  energy  into  the  electron  gas  becomes  very 

® —6 

complex.  For  f^ > 10  , however,  the  deposition  can  again  be  represented 

simply  by  (see  Appendix  A) : 


oc  , 


(4-12) 


This  aleration  of  form  is  of  more  than  passing  interest  since  fractional 
ionizations  of  10  and  above  occur  quite  commonly  at  altitudes  i 200  km. 


Using  Equation  4-12,  we  can  now  write  the  heating  rate  per 


electron  as 


d,?  /dt 

An  ' ^ ^ 

e 


(4-13) 
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Equation  4-13  tells  us  that  while  the  total  energy  deposition  rate  in  the 
electron  gas  increases  as  the  electron  density  increases,  the  heating  rate 
per  electron  actually  decreases  as  the  electron  density  increases.  This  is 
in  contrast  to  the  previous  result  (Equation  4-10)  where  it  remained 
constant.  Similarly,  the  steady  state  temperature,  will  be  a function 

of  n , increasing  with  decreasing  n . 

C 0 

The  dynamic  solution  to  Equation  4-1  for  the  electron  temperature, 
under  conditions  where  Equation  4-12  is  valid,  is  represented  schematically 
in  Figure  4-2.  We  see  that,  as  in  Figure  4-1,  there  is  an  initial  tempera- 
ture rise  but,  unlike  that  in  Figure  4-1,  the  temperature  ultimately  reaches 
a maximum  and  then  subsequently  decays.  This  behavior  can  be  understood 
qualitatively  by  writing  the  solution  to  Equation  4-1  as 


T (T  -T  )e  + T 
e e ss  ss 

o 


(4-14) 


Quasi -Steady  State  Temperature  T 


•Electron  Temperature 


A 


Neutral  Gas  Temperat 


Figure  4-2.  Behavior  of  Tg  and  Tgj  using  Equation  4-12  and  4-13. 
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where  x is  determined  from  Equation  4-11  and  T is  the  solution  to  | 

ss  j 

H = C(Tgg)  . (4-15)  S 

j 

A schematic  plot  of  is  also  shown  in  Figure  4-2.  j 

I 

Initially,  in  Figure  4-2,  the  electron  temperature  rises  because  | 

the  bombarding  flux  heats  the  gas.  This  is  represented  by  the  exponential  I 

-4  1 

term  in  Equation  4-14.  The  rise  time,  x,  can  vary  from  10  sec  at  100  km 
to  seconds  at  200  km.  During  this  rise  time,  however,  the  electron  density 
is  increasing.  This  increase  results  in  a decrease  of  the  steady  state 
electron  temperature  and,  subsequently,  once  T reaches  T , to  a decay 
of  T^.  The  decay  of  T^  is  almost  entirely  governed  by  the  time  constant 
for  electron  density  buildup  which  is  usually  longer  than  the  time  constant 
of  Equation  4-11. 

Figures  4-3  and  4-4  present  the  results  of  calculations  at  severjJ 
different  altitudes  (using  Equation  A-28  of  Appendix  A for  the  partition 
function)  for  two  different  input  'spectra.  Figure  4-3  shows  the  electron 
temperature  at  100,  132  and  194  km,  as  functiohs  of  time,  under  electron 
flux  conditions  corresponding  to  the  '75  Multi  arc  (Section  2).  We  see 
that  the  qualitative  behavior  expected  (Figure  4-2)  is  indeed  borne  out. 

Further,  we  see  that  at  quite  low  altitudes  (132  km)  the  transient  electron 
temperature  can  easily  reach  large  values  (2000°K)  before  decaying  away. 

Figure  4-4  presents  the  same  information  for  the  background  drizzle  appro- 
priate to  the  '74  Multi  event  (Reference  4).  The  energy  input  from  this 
spectrum  is  much  less  than  for  the  '75  Multi  arc,  but  even  here  we  see  ■ 

some  transient  behavior  as  in  Figure  4-2. 


HEATING  BY  N(^D)  QUENCHING 

As  mentioned  earlier,  deexcitation  of  excited  atomic  states  by  the 
electron  gas  will  also  contribute  to  the  heating  rate.  The  reactions  are 
typically  of  the  type 

e + X*  ^ X + e*  (4-16) 

where  the  recoil  electron  can  be  approximated  as  having  an  energy  equal  to 
the  energy  of  the  excited  state  x*.  The  additional  energy  carried  by  the 
"hot"  electron  e*  will  now  be  deposited  in  the  gas.  Some  of  this  energy 
will  go  into  thermal  excitation  of  the  electron  gas.  The  partition  of  the 
energy  deposited  by  these  low  energy  electrons,  incorporated  in  the  ARCTIC 
code,  is  given  in  detail  in  Appendix  B, 

For  most  states,  the  additional  heating  due  to  reactions  like 

2 

4-16  is  very  small  and  can  be  neglected.  However,  for  N(  D)  the  de- 
excitation reaction 

N(^D)  + e-J^N(''s)  + e(2.38  eV)  (4-17) 

can  play  a significant  role  at  high  altitudes.  In  fact,  in  Appendix  B we 
show  that  virtually  all  of  the  energy  of  the  hot  electron  goes  into  thermal 
energy.  Therefore,  the  rate  at  which  energy  is  being  added  to  the  thermal 
electron  gas  is 

(d^/dt)j^^2D)  ^ 2-38K[N(2D)]ng  (eV  cm‘^  sec'^)  (4-18) 

2 

and  the  heating  rate  due  to  deexcitation  of  N(  D)  can  be  written,  using 
Equation  4-7,  as 


■1 
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In  Figure  4-5  we  have  plotted  the  fraction, 

2 

heating  rate  that  arises  from  deexcitation  of  N(  D) . 


f,  of  the  total 
It  is  defined  as 


H(N(^D)) 

H(N^D))  + H (bombarding  flux) 


(4-20) 


The  results  in  Figure  4-5  are  for  the  '75  Multi  arc  at  an  altitude  of  194  km. 

They  show  that  f varies  from  insignificant  amounts  at  times  less  than 

10  sec  to  a value  in  excess  of  20  percent  of  the  total  heating  rate  in 

steady  state.  We  would  expect  from  Equation  4-20  that  it  would  follow  the 
2 

N(  D)  concentration  closely,  and  indeed  it  does,  as  can  also  be  seen  in 
Figure  4-5.  Of  course  the  bombarding  flux  is  still  the  ultimate  energy  source 
for  electron  heating.  What  we  have  shown,  however,  is  that  indirect 
processes  can,  at  times,  become  significant  in  determining  the  electron 
temperature . 

COMPARISONS  WITH  SRI  RADAR  DATA 


The  foregoing  procedures  were  incorporated  into  the  ARCTIC  code. 
New  runs  were  then  made  for  a variety  of  auroral  input  conditions.  The 
resulting  electron  temperatures,  as  functions  of  altitude,  are  shown  by  the 
solid  curves  in  Figure  4-6.  Also  shown  (dashed  curve)  are  the  neutral  gas 
temperature  and  several  data  points  inferred  from  the  SRI  incoherent  back- 
scatter  radar.  The  only  common  auroral  events  that  can  be  used  for  compari- 
son are  those  on  27  March  1973  ('73  Multi)  and  25  February  1974  ('74  Multi). 

For  the  '73  Multi  event,  the  calculations  give,  for  an  altitude 
of  175  km,  arc  and  drizzle  temperatures  of  1360  "K  and  900  °K,  respectively, 
compared  with  a radar  value  of  1250  ®K.  For  the  '74  Multi  event,  the  calcu- 
lations give,  for  an  altitude  of  170  km,  arc  and  drizzle  temperatures  of 
1600  ®K  and  1030  ®K,  respectively,  compared  with  a radar  value  of  1300  ®K. 


d 


58 


Time  (sec) 


Figure  4-5 


. Fraction,  f,  of  total  electron  heating  rate  at  194  km  due  to 
deexcitation  of  NC^d)  after  onset  of  '75  Multi  arc  flux. 


[x]  (Uncertain  by  SOX) 
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25  ergs  cm  sec  ; <E>  » 1.2  kev 


Tg(2)  - Calculated  (ARCTIC);  1973  Multi-arc;  141  ergs  cm'^sec"' 

Tg(3)  - Calculated  (ARCTIC);  1974  Multi-arc;  285  ergs  cm"^sec'^ 

Tg(4)  - Calculated  (ARCTIC);  1975  Multi-arc;  59  ergs  cm'^sec’’ 

Tj(5)  - Calculated  (ARCTIC);  1974  Multi-drizzle;  3.9  ergs  cm'^sec"' 
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Figure  4-6.  Electron  temperature  calculations  compared  with  SRI  radar  data. 
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Since  the  radar  measurements  may  correspond  to  conditions  either 
inside  or  outside  the  arc,  or  to  a combination  of  both,  it  is  perhaps 
reasonable  that  the  data  should  fall  between  the  calculated  arc  and  drizzle 
values,  as  indeed  it  does  for  both  of  these  events. 
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SECTION  5 

NO  IN  THE  AURORAL  ENVIRONMENT 


PRELIMINARY 

Enhancement  of  the  NO  concentration  in  an  auroral  environment  has 
been  reported  by  a number  of  workers  (see,  for  example.  References  18,  19,  and 
20) . In  an  earlier  report  (Reference  2),  we  also  inferred  enhanced  values  for 
[NO]  based  on  ion  mass  spectrometer  measurements  for  the  ICECAP  event  of 
27  March  1973.  Furthermore,  it  is  difficult  to  account  for  the  ICECAP  data 
near  5.3  pm  unless  we  attribute  it  to  the  presence  of  an  enhanced  NO  profile. 

With  the  exception  of  the  neutral  mass  spectrometer  measurements 
of  Zipf,  et  al.  (Reference  18) , most  measurements  on  NO  concentrations  in 
auroras  are  indirect.  We  have  taken  the  relevant  ICECAP  data  and  have 
attempted  to  infer  NO  profiles  on  the  basis  of  two  different  methods,  neither 
of  which  is  very  reliable.  However,  although  the  uncertainties  are  large, 
the  results  unanamously  support  the  presence  of  NO  concentrations  consider- 
ably larger  than  those  reported  at  lower  latitudes.  The  results  of  this 
work  are  reported  here. 

[NO]  INFERRED  FROM  SLOPE  OF  5,3-ym  ZENITH  INTENSITY  PROFILES 

The  basis  for  this  method  of  inferring  [NO]  profiles  is  the  observa- 
tion that  most  of  the  ICECAP  data  on  zenith  intensities  at  5.3  ym  are,  for  a 
given  event,  essentially  independent  of  whether  the  sight  path  intersects  an 
auroral  arc  or  not.  An  exception  to  this  rule  may  occur  only  for  very 
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intense  events.  The  only  way  we  have  been  able  to  account  for  these  observa- 
tions is  to  assume  that  the  NO  concentration  is  uniformly  enhanced  both 
outside  and  inside  the  arc.  Excitation  of  the  NO(v=l)  state  then  occurs 
by  a combination  of  the  atom  interchange  reaction 


N0(v=0)  + 0' 


-►  NO  (v=l)  + 0 


(5-1) 


and  earthshine  scatter  by  NO.  For  sufficiently  enhanced  concentrations  of 
NO,  this  emission  dominates  the  NO  chemiluminescence  contribution,  except 
possibly  in  very  intense  arcs. 


To  the  extent  that  this  view  is  correct,  one  can  then  infer  an 
NO  concentration  by  noting  that  it  is  proportional  to  the  volume  emission 
rate  in  the  NO  fundamental  band  near  5.3  ym,  especially  outside  arcs. 
Mathematically,  since  ij)(z),  the  volume  emission  rate  at  altitude  z,  is 
related  to  I(z),  the  zenith  intensity,  by 

00 


1(2) 


^dz 


(5-2) 


we  have 


<(>(z)  = - 


dl(z) 


dz 


(5-3) 


But  under  conditions  where  the  excitation  of  N0(v=l)  occurs  entirely  by 
atom  interchange  + earthshine  scatter,  we  have 


(5-4) 


<p  = (k[0]  +n)[N0] 

where  k is  the  rate  coefficient  for  Reaction  5-1  and  n is  the  earthshine 
excitation  rate  per  NO  molecule.  Equations  5-3  and  5-4  then  lead  to  the 


result 


fNOl  = ~ (dl/dz) 

^ ^ (k[0]  + n) 


(5-5) 
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Normally,  the  measuring  instrument  has  limited  bandwidth  so  that 
only  a fraction,  f,  of  the  emitted  photons  are  detectable.  Thus,  if  (dl/dz)' 
is  the  slope  of  the  measured  intensity  profile.  Equation  5-5  can  be  written 
as 


[NO] 


-(dl/dz)  " 

f(k[o] 


(5-6) 


Equation  5-6  has  been  used  to  estimate  NO  profiles  for  three 
different  auroral  events.  The  results  are  shown  in  Figure  5-1.  They  are 
based  on  radiometer  data  obtained  on  descent  through  drizzle  regions  by 
Rockets  A18. 205-1  ('73  Multi)  and  IC519.07-1B  ('75  Multi)  and  by  Rocket 
IC507.11-2A  during  the  intense  auroral  breakup  event  of  10  March  1975. 
Values  used  for  the  parameters  appearing  in  Equation  5-6  are  (see  Reference 
4):  f = 0.75,  0 = 3.38X10"^,  k = 5 x lo'^^  .j,0.5g-3196/T^ 

Figure  5-1  (dashed  curve)  is  a mid  latitui,  profile  obtained  from  Refer- 
ences 21  and  22. 


The  uncertainty  associated  with  the  profiles  shown  in  Figure  5-1 
is  large  because  of  the  difficulty  in  accurately  determining  the  slope, 
(dl/dz)',  of  the  measured  intensity  profiles,  and  because  of  uncertainties 
in  the  values  of  k and  q.  Factors  of  2 or  more  in  uncertainty  of  the 
inferred  NO  values  would  not  be  unreasonable. 

Application  of  this  method  to  the  drizzle  environments  appropriate 
to  the  descent  airspace  of  the  '73  and  '75  Multi  rockets  seems  reasonable 
since  NO  chemiluminescence  should  be  small.  However,  one  might  question 
its  use  for  the  intense  auroral  breakup  event  of  10  March  1975  in  which 
chemiluminescence  would  be  expected  to  contribute  significantly.  The  only 
justification  is  based  on  the  following  facts  (see  Section  3);  (1)  the 
measured  zenith  intensity  in  the  radiometer  band,  centered  near  5.6  ym, 


64 


Alt 


1 n 


Figure  5-1.  Auroral  NO  concentrations  inferred  from  slope  of 
5.3-ym  radiometer  data. 
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is  nearly  the  same  on  rocket  ascent  and  descent,  (2)  the  computed  chemi- 
luminescence is  too  small  to  account  for  the  data,  and  (3)  the  CVF  spectra 
support  the  view  that  the  dominant  emission  is  the  (1,0)  band  of  NO. 

[NOl  INFERRED  FROM  ION  MASS  SPECTROMETER  MEASUREMENTS 

It  can  be  shown  (see,  for  example.  Reference  19)  that  in  steady  state 
the  NO  concentration  can  be  determined  from  measured  ratios  of  the  ion 
concentrations  by  the  expression 

M • r 

^2  ®\[02  V 2 \[02  1/  ’^2  Mt02  1/  2 

Here,  the  k's  are  rate  constants  for  the  reactions: 


N2‘^  + 0 - 

^1 

-—^N  + 

NO* 

; kj  = 

-9  -0  44 

1.7  X 10  ^ T ■ ^ 

02"’  + NO 

^2 

--^^2 

+ NO* 

» ^^2  “ 

6.3  X 10"^° 

NO*  + e - 

-^-►N  + 

0 

; k3  = 

3.5  X 10"^  T 

e 

0*  + N2  - 

4 't- 

— ^NO 

+ N 

3.6  X 10’^^  T"^ 

N*  . O2  - 

k 

— ^NO* 

+ 0 

' h = 

3 X 10’^® 

and  n^  is  the  electron  density. 


* Actually,  we  don't  know  if  the  particle  precipitation  on  rocket  descent 
was  significantly  less  than  on  ascent  because  the  SRI  radar  was  not 
operated  in  a meridian  scan  mode  during  the  rocket  flight  but,  rather, 
was  pointed  at  the  anticipated  100-km  altitude  rocket  entry  point. 
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We  have  used  Equation  5-7,  together  with  measured  (Reference  23) 
ion  mass  spectrometer  values  for  the  ion  concentration  ratios,  and  retarding 
potential  analyzer  values  for  n^,  to  determine  [NO]  profiles  on  ascent  and 
descent  of  the  '73  Multi  Rocket  A18. 205-1  on  27  March  1973.  Values  for  the 
neutral  species  concentrations  were  taken  from  the  1965  Cira  mean  model 
(Model  5,  hour  0 above  120  km),  and  values  for  the  electron  temperature, 
obtained  from  the  ARCTIC  code,  are  shown  by  the  appropriate  curves  in 
Figure  4-6.  The  results  are  shown  in  Figure  5-2. 

The  jagged  appearance  of  the  curves  in  Figure  5-2  arises  from  the 
fact  that  in  Equation  5-7  we  used  the  actual  values  of  the  ion  ratios  that 
exhibit  jagged  profiles  due  to  spin  modulation  of  the  rocket. 

A similar  calculation  has  been  performed  by  Narcisi  and  Swider 
(Reference  20)  for  the  ascent  portion  of  the  rocket  flight,  using  smoothed 
versions  of  their  ion  concentration  ratios.  Their  results,  for  the  Cira 
1965  atamosphere,  match  ours  at  110  km,  but  fall  off  faster  with  increasing 
altitude.  We  attribute  the  difference  to  their  use  of  different  rate 
constants,  to  their  equating  the  electron  and  neutral  temperatures  and, 
perhaps,  to  their  use  of  somewhat  different  neutral  species  concentrations. 

Uncertainty  in  the  calculated  NO  profile,  using  Equation  5-7, 
can  arise  for  at  least  two  reasons.  First,  when  the  terms  in  Equation  5-7 
are  of  comparable  magnitude,  especially  the  first  two,  small  variations  in 
the  magnitude  of  their  factors  produce  large  changes  in  the  resulting  NO 
concentration.  Second,  the  equation  is  valid  only  in  the  steady  state,  a 
condition  that  does  not  generally  prevail  in  a dynamic  auroral  environment. 

Also  shown  in  Figure  5-2  (smooth  solid  curve)  for  comparison,  is 
the  descent  NO  profile  from  Figure  5-1  obtained  by  taking  the  slope  of 
the  5.3-ym  radiometer  data.  Over  the  altitude  range  from  105  to  140  km. 
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NO  Concentration  (cm"'^) 

Figure  5-2.  NO  Concentration  inferred  from  auroral  data  of  27  March  1973 
(Rocket  A1 8. 205-1). 
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^ the  latter  profile  is  seen  to  be  smaller  than  the  one  obtained  by  the  ion 

ratio  method  by  factors  ranging  from  about  10  at  110  km  to  about  3 or  so  at 
140  km. 

i. 

Two  mid  latitude  profiles  (widely  different)  are  also  shown  in 
' Figure  5-2.  In  all  cases,  the  auroral ly- inferred  profiles  are  considerably 

enhanced  over  those  at  mid  latitudes. 

SUMMARY 

] 

Although  fraught  with  considerable  uncertainty,  two  different 
methods  have  been  used  to  infer  [NO]  profiles  from  ICECAP  data.  The 
results  are  all  supportive  of  considerable  auroral  enhancement  of  [NO], 
especially  at  altitudes  in  the  100-  to  120-km  region. 
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SECTION  6 

ICECAP  SUfflARY  AND  CONCLUSIONS 


i 


PRELIMINARY 

This  section  summarizes  the  principal  results  and  conclusions  of 
our  work,  as  embodied  in  this  and  earlier  reports  (References  1 through  4) , 
related  to  the  interpretation  and  understanding  of  ICECAP  data.  Emphasis 
has  been  on  those  auroral  events  for  which  coordinated  measurements  were 
taken  with  a variety  of  instruments,  especially  the  short  wavelength  infrared 
(SWIR)  detectors.  Not  all  aspects  of  the  ICECAP  data  have  been  studied,  and 
no  attempt  is  made  here  to  summarize  results  from  the  entire  ICECAP  program, 
or  conclusions  reached  by  other  workers. 

One  of  the  main  purposes  of  ICECAP  was  to  measure  SWIR  radiation 
from  an  aurorally  excited  atmosphere  and  thereby  test  the  validity  of  models 
used  in  computer  codes  to  predict  corresponding  emissions  in  a nuclear  en- 
vironment. In  this  sense  the  program  was  at  least  partially  successful. 
Complete  success  was  hampered  by  a number  of  factors,  not  the  least  of  which 
is  the  fact  that  the  experiments  had  to  be  carried  out  in  the  uncontrolled, 
dynamic,  auroral  environment.  It  is  not  possible,  in  such  experiments,  to 
measure  the  input  particle  flux  spectriim  over  a large  enough  spatial  region 
and  for  sufficiently  long  periods  of  time  to  specify  adequately  the  input 
conditions  necessary  for  accurate  calculation  of  the  resulting  optical/IR 
radiances.  Ambiguities  related  to  excitation  mechanisms  and  radiating 
species  also  exist  because  of  instrumental  limitations,  especially  with 
respect  to  spectral  coverage  and  wavelength  resolution.  In  some  cases, 
radiance  data  at  a given  wavelength,  measured  by  two  different  instruments. 
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may  be  sufficiently  different  to  render  the  field  values  very  uncertain.  So, 
for  these  and  other  reasons,  uncertainty  still  remains  with  respect  to  in- 
terpretation of  some  of  the  data,  and  complete  verification  of  the  SWIR 
models  is  not  yet  possible. 

Our  analysis  and  interpretation  of  the  data  was  greatly  facilitated 
by  the  development  of  the  ARCTIC  code  that  inputs  the  measured  spectral  flux 
of  auroral  electrons  and  computes  in  detail  the  energy  deposition  and  par- 
tition as  well  as  the  subsequent  atmospheric  chemistry  and  optical/IR 
radiation.  Many,  although  not  all,  of  our  conclusions  are  based  on  compari- 
sons between  these  calculations  and  the  field  data. 

PRINCIPAL  RESULTS 

The  following  summarizes  our  results  according  to  the  different 
topics  covered.  To  simplify  reference  to  the  different  auroral  events,  we 
have  classified  them  numerically  according  to  the  scheme  shown  in  Table  6-1. 

Table  6-1.  Classification  of  auroral  events  studied. 


Event  Date  and  Description 

Rocket 

27  March  1973  (discrete  arc) 

24  March  1973  (auroral  breakup) 

25  February  1974  (discrete  arcs) 

12  March  1974  (discrete  arcs) 

10  March  1975  (auroral  breakup) 

A18. 205-1:  Black  Brant  ( '73  Multi) 

AID. 205-2;  Paiute  Tomahawk 

A18. 219-1;  Black  Brant  ( '74  Multi) 

IC519.07-1B;  Sergeant  Hydac  ('75 Multi) 

IC507.11-2A:  Nike  Hydac 

iVRlVIiK'^vi 


5.3-ym  Radiation 

Radiometer  and/or  CVF  spectrometer  data  from  each  of  the  five 
events  shown  in  Table  6-1  were  analyzed  and  compared  with  calculations.  The 
results  can  be  enumerated  as  follows: 

(1)  NO  chemiluminescence  cannot  account  for  any  of  the  data.  Its 
magnitude  (too  small) , its  sensitivity  to  rapid  changes  in  bombardment  con- 
ditions, and  its  spectral  distribution,  are  not  in  accord  with  the  observa- 
tions. 

(2)  Although  generally  differing  in  magnitude  from  event  to  event,  the 
zenith  intensity  profile  near  5.3  ym,  for  any  given  event,  is  nearly  the 
same  on  rocket  ascent  and  descent,  even  though  the  instantaneous  bombard- 
ment conditions  may  be  quite  different  in  the  upleg  and  downleg  regions.  A 
good  example  of  this  behavior  is  shown  in  Figure  6-1  for  auroral  Event  2 

O 

(Table  6-1) . Although  the  zenith  3914A  intensity  on  ascent  saturated  at 
200  kR,  the  intensity  on  descent  did  not  exceed  about  40  kR  (see  Reference  4, 
Figure  3-2).  However,  at  least  above  95  km  where  we  have  descent  data,  the 
5.3-ym  intensity  is  virtually  the  same  on  ascent  and  descent. 

(3)  The  limited-resolution  spectral  data  are  reasonably  consistent 
with  (1,0)  band  radiation  from  NO. 

(4)  Excitatirn  of  N0(v  = 1)  by  a combination  of  the  atom  interchange 
reaction  N0(v  = 0)  +0'^N0(v  = l)  +0  and  earthshine  scatter  appears  to  be 
capable  of  accounting  for  most  of  the  data.  However,  these  mechanisms  do 
require  the  presence  of  considerably  enhanced  [NO]  profiles  that  differ 
from  event  to  event.  Complete  verification  of  this  model  would  entail  an 
independent  measurement  of  the  [NO]  profiles  as  well  as  a measurement  of  the 
rate  constant  for  the  atom  interchange  reaction. 

The  data  at  5.3  ym  have  produced  one  of  the  more  surprising  out- 
comes from  the  ICECAP  program.  The  results  have  implications  for  modeling 
in  the  nuclear  environment  that  will  be  mentioned  later. 
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Peak  Spectral  Radiance  (MRytn" 


Rocket  Altitude  (km) 


Figure  6-1.  Zenith  peak  spectral  radiance  near  5.3  ym 
(Paiute  Tomahawk  A10-205-2). 
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2.7-vtm  Radiation 


I 

js' 


i 

I 


f 


i 

! 

I 


i 


t 

' I 

: ' i 


The  only  usable  rocket  data  at  2.7  pm  from  ICECAP  that  we  are 
aware  of  were  taken  by  CVF  instruments  during  Events  2 and  3 (Table  6-1) . 

The  CVF  data  for  Event  4 are  below  the  instrument  noise  level.  From  study 
of  this  very  limited  amount  of  data,  and  comparison  with  calculations,  we 
find  the  following; 

(1)  For  Event  2,  comparison  between  the  measured  zenith  intensity  at 

O 

3914A  and  that  near  2.7  ym,  leads  to  a fluorescence  efficiency  for  the  2.7-ym 
band  of  about  0.5  percent.  This  is  consistent  with  analyses  of  aircraft 
data  reported  by  Kofsky  (Reference  25) . It  is  also  consistent  with  our  code 
calculations  of  NO  chemiluminescence  that  provides  for  efficiencies  at  2.7  ym 
between  about  0.5  and  1 percent. 

(2)  For  Event  2,  the  magnitude  of  the  measured  intensity  above  95-km 

altitude,  and  the  shape  of  the  spectral  feature  near  2.7  ym,  are  reasonably 

consistent  (to  within  a factor  of  2 on  intensity)  with  our  model  of  NO 

2 

chemiluminescence  that  provides  for  about  one  2.7-ym  photon  per  N(  D)  +0^ 
reaction. 

(3)  For  Event  3,  NO  chemiluminescence  is  insufficient  to  account  for 
the  magnitude  of  the  2.7-ym  data  available  to  us.  In  fact,  calculation  of 
the  fluorescence  efficiency  at  2.7  ym,  obtained  by  comparing  the  CVF  data 

O 

with  the  3914A  data,  leads  to  values  as  high  as  50  percent.  Therefore,  it 
would  appear  that  either  there  is  a calibration  problem  with  the  CVF  instru- 
ment on  this  flight,  or  else  we  don't  understand  the  radiation  mechanism. 

For  Event  3,  the  data  referred  to  above  were  taken  on  rocket 
descent  about  one  minute  after  passage  of  a strong  arc.  It  has  been 
suggested  (Reference  26)  that  a long-life  radiator  (perhaps  CO2)  is  operating 
near  2.7  ym  in  the  auroral  environment.  If  so,  this  could  account  for  the 
observed  magnitude  of  the  emission  which  would  have  been  excited  by  the 
earlier  passage  of  the  arc.  However,  the  relatively  slow  decay  of  the 
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intensity  with  increasing  altitude  above  about  90  km,  compared  with  that  at 
4.3  ym,  does  not  favor  CO2  as  the  likely  radiator. 

In  summary  then,  except  for  Event  3,  the  weight  of  the  very  limited 
data  (aircraft  + one  rocket  flight)  suggest  that  NO  chemiluminescence  can 
account  reasonably  well  for  the  ICECAP  observations  near  2.7  ym. 

4.3-ym  Radiation 

Radiometer  and/or  CVF  spectrometer  data  for  the  first  four  events 
shown  in  Table  6-1  were  studied  and  compared  with  calculations  involving 
CO2  vibraluminescence  and  NO*  chemiluminescence.  The  results  are  as  follows: 

(1)  The  data  are  all  reasonably  consistent  with  CO2  vibraluminescence 
as  the  main  source  of  4.3-ym  radiation  in  the  auroral  environment.  Both  the 
magnitude  of  the  intensity  and  the  shape  of  the  spectral  features  near  4.3  ym 
support  this  conclusion.  The  emission  is  consistent  with  the  presence  of  a 
long-life  radiator.  Removal  of  an  intense  bombardment  source  (arc)  does 

not  immediately  reduce  the  emission  intensity  as  would  be  the  case  if  the 
source  were  primarily  NO*  chemiluminescence.  For  this  reason,  however, 
detailed  agreement  between  calculations  and  the  field  data  is  difficult  to 
achieve  because  transport  effects,  due  to  winds  and  diffusion,  have  time  to 
modify  the  spatial  distribution  of  the  vibrational ly  excited  nitrogen  and, 
hence,  the  CO2  volume  emission  rate.  Knowledge  of  the  bombardment  history 
(for  minutes  to  hours,  depending  on  the  altitude)  is  also  essential  for 
correct  determination  of  the  vibraluminescent  intensity. 

(2)  Unquenched  NO*  chemiluminescent  emission  near  4.3  ym  is  generally 
in  poor  agreement  with  the  data,  both  with  respect  to  the  magnitude  of  the 
intensity  and  to  the  shape  of  the  spectral  features.  Also,  as  mentioned 
above,  the  observed  radiation  persists  after  removal  of  the  bombardment 
source,  and  this  is  not  consistent  with  NO*  chemiluminescence. 
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Data  from  Event  4 has  been  used  to  show  that  a rate  constant  as 
-93-1  + 

large  as  1 x 10  cm  sec  for  quenching  NO  emission  by  N2  is  required. 

The  likely  mechanism  is  v - v transfer  from  NO  to  N2  for  which  there  is 
very  close  energy  resonance. 

NO  Concentration 

Values  for  the  NO  concentration,  enhanced  by  up  to  2 orders  of 
magnitude  over  mid  latitude  values  near  105  km,  have  been  inferred  from 
Events  1,  4,  and  5 (Table  6-1).  Even  larger  values  for  Event  2 are  implied 
by  the  data  in  Figure  6-1. 

Some  of  these  inferred  values  (for  Event  1)  are  based  on  ion  mass 
spectrometer  data.  The  others  are  based  on  the  assianption  that  most  of  the 
5.3-ym  radiation  is  excited  by  the  atom  interchange  reaction + earthshine 
scatter  from  NO.  To  the  extent  that  this  is  true,  the  data  suggest  that  the 
altitude  profile  of  NO  is  fairly  uniform  over  a large  spatial  region  (perhaps 
a few  hundred  km)  surrounding  a given  auroral  form.  The  degree  of  enhance- 
ment s^pears  to  be  related  to  the  cumulated  energy  deposited  over  some  period 
of  time,  although  there  does  seem  to  be  a correlation  with  the  intensity  of 
the  specific  event  as  well.  To  our  knowledge,  no  one  has  attempted  to  take 
detailed  account  of  the  bombardment  history,  for  a given  event,  over  a period 
of  perhaps  hours  or  days,  and  perform  a chemistry  calculation  to  see  if  the 
required  NO  values  can  indeed  be  attained.  Part  of  the  difficulty  in  making 
such  a calculation,  however,  is  the  unknown  influence  of  transport  effects 
that  have  ample  time  to  operate  and  change  the  NO  distribution. 

IMPLICATIONS  TO  THE  NUCLEAR  ENVIRONMENT 

Application  of  the  foregoing  results  to  a nuclear  environment  are 
fairly  clear,  at  least  for  the  5.3-  and  4.3-pm  bands. 


For  the  5.3-yin  band,  the  models  should  add  the  contribution  from 
NO  radiation  due  to  excitation  by  the  atom  interchange  reaction  and  by  earth- 
shine  scatter.  We  have  not  studied  the  effect  of  this  addition  in  detail, 
but  have  made  estimates  for  two  situations  involving  high-altitude  nuclear 
detonations.  In  one  case,  for  a sight  path  along  which  there  was  no  signi- 
ficant temperature  increase,  the  effect  was  to  approximately  double  the 
calculated  radiance  due  to  NO  chemiluminescence  alone.  In  the  other  case, 
for  sight  paths  external  to  the  fireball,  but  along  which  there  was  signi- 
ficant heating,  the  atom  interchange  contribution  was  greater  by  up  to  2 
orders  of  magnitude  than  the  combined  effects  of  NO  thermal  and  chemi- 
luminescence radiation. 

For  the  4.3-ym  band,  CO^  vibraluminescence  seems  clearly  established 
as  the  dominant  mechanism,  at  least  in  regions  not  heated  sufficiently  for 
thermal  emission  to  dominate.  However,  modeling  for  the  nuclear  environ- 
ment is  complicated  in  two  ways.  First,  as  interpretation  of  the  ICECAP 
data  by  Kumer  (Reference  13)  has  shown,  trapping  of  the  4.3-ym  photons  at 
altitudes  below  about  100  km  is  important.  This  is  especially  significant 
in  the  nuclear  environment  because  much  of  the  X-ray  energy  is  deposited  at 
altitudes  near  80  km.  Second,  because  of  the  long  lifetime  of  the  radiation 
process,  the  effects  due  to  transport  of  the  vibrationally  excited  nitrogen 
by  winds  and  diffusion  cannot  be  ignored.  The  models  must  make  some  pro- 
vision for  incorporating  these  effects  if  reasonably  reliable  predictions 
are  to  be  made.  This  is  particularly  so  for  the  case  of  a multiburst 
environment . 

For  the  2.7-ym  band,  the  results  are  still  in  some  doubt.  At  the 
present  time,  we  would  recommend  retention  of  a chemiluminescent  model  in 
which  about  0.5  percent  of  the  deposited  energy  is  radiated  in  the  first 
overtone  band  of  NO.  Further  work  is  needed  to  clarify  whether  or  not 
some  other  molecule,  or  excitation  mechanism,  is  contributing  significantly 
to  this  band,  at  least  at  altitudes  above  about  90  km. 
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SECTION  7 

COMMUNICATIONS>RELATED  WORK 


INTRODUCTION 


With  DNA's  increased  emphasis  on  field  programs  related  to  satellite 
communications,  especially  the  Wideband  Satellite  experiment,  a small  part  of 
our  present  effort  was  devoted  to  work  in  this  area. 

Satellite  communication  problems  arise  largely  from  radiowave 
scintillation  caused  by  striations  or  other  irregularities  in  the  F-region 
ionization.  Successful  prediction  of  the  irregularity  structure  and  scin- 
tillation effects  in  a nuclear  environment  depends  on  a good  theoretical 
understanding  of  how  the  irregularities  develop.  Hopefully,  a major  step 
toward  this  understanding,  at  least  in  the  non  nuclear  case,  will  result 
from  the  proper  interpretation  of  the  Wideband  Satellite  transmissions 
through  both  the  equatorial  and  high  latitude  scintillation  zones. 


Two  aspects  of  the  problem,  especially  as  it  bears  on  the  inter- 
pretation of  the  Wideband  data,  were  addressed.  One  relates  to  the  question 
of  how  important,  in  the  auroral  environment,  the  perturbing  effects  due  to 
structure  in  particle  precipitation  may  be  in  masking  the  instability- 
produced  inhomogeneities  in  the  F layer.  The  other  relates  to  the  question 
of  how  to  measure  the  F-region  neutral  wind  component  normal  to  the  earth's 
magnetic  field  at  equatorial  latitudes.  This  wind  component  is  the  driving 
force  in  the  so-called  "ion-neutral  slip"  (or  "gradient  drift")  instability 
mechanism  that  is  believed  to  generate  the  striation  phenomenon.  Although 
our  efforts  in  these  two  areas  were  not  very  fruitful,  they  are  briefly 
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AURORAL  IONIZATION  EFFECTS  DUE  TO  STRUCTURED  PARTICLE  FLUXES  i 

Interpretation  of  Wideband  Satellite  scintillation  data  from  the  j 

auroral  region,  in  terms  of  instability-generated  F-layer  inhomogeneities,  | 

could  be  rendered  erroneous  if  effects  due  to  structure  (temporal  or  spatial) 
in  the  auroral  particle  flux,  were  to  dominate.  In  this  context,  we  first 
note  that  direct  effects  on  F-layer  ionization  from  a structured  particle 
flux  should  be  insignificant  because  of  the  long  response  time  of  the 
ionosphere  above  250  km.  In  the  E layer,  however,  ionization  irregularities 
would  not  be  expected  to  wash  out  because  of  the  much  shorter  response  time. 

We  are  not  referring  here  to  structure  associated  with  discrete 
arcs.  For,  although  auroral  arcs  may  produce  severe  radiowave  scintillation  | 

as  the  propagation  path  scans  across  them,  their  presence  in  the  Wideband 
experiment  would  surely  be  recognized,  either  from  coordinated  optical 
measurements,  or  from  their  characteristic  scintillation  signature.  Rather, 
the  concern  was  over  the  possibility  of  frequent  occurrence  of  random,  large 
amplitude,  fluctuations  of  the  incoming  particle  flux.  Such  occurrences 
would  then  be  reflected  in  corresponding  fluctuations  in  the  E-layer  ioniza-  j 

tion  that  might  be  misinterpreted  in  terms  of  F-layer  structure.  ^ 

We  conducted  a fairly  exhaustive  literature  survey  to  uncover 
what  is  known  about  such  random  flux  variations.  The  outcome  of  that  survey 
shows  that,  for  the  most  part,  the  particle  precipitation  in  regions  away 
from  discrete  auroral  forms  is  quite  uniform.  The  presence  of  large  amplitude  i 

landom  fluctuations  is  rare.  Consequently,  we  conclude  that  the  probability 
that  structured  particle  fluxes  will  interfere  with  auroral  Wideband  data 

interpretation  is  small.  ? 

.1 
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ELECTRON  BEAM  EXCITATION  OF  F-REGION  ATOMIC  OXYGEN 

As  was  mentioned  in  the  Introduction,  the  force  driving  the  in- 
stability mechanism  believed  responsible  for  F-region  striations  is  the 
neutral  wind  component  normal  to  the  earth's  magnetic  field.  A suggested 
method  for  measuring  this  wind  component  is  to  excite  one  or  more  lines  in 
the  visible  spectrum  of  atomic  oxygen  by  an  in  situ  electron  beam  and  to 
measure  the  doppler  shift  in  the  lines  by  utilizing  a ground-based  Fabry 
Perot  interferometer.  The  following  paragraphs  describe  the  results  of  a 
study  made,  using  the  ARCTIC  code,  to  assess  the  utility  of  such  an  experi- 
ment at  equatorial  latitudes.  They  show  why  we  feel  that  the  technique 
is  not  feasible. 

Deposition  and  Partition  by  325  eV  Electrons  at  300-km  Altitude 

The  ARCTIC  code  was  run  to  determine  the  deposition  volume  and 
the  partition  of  the  energy  for  an  electron  beam  of  various  initial  energies 
at  an  altitude  of  300  km  at  the  magnetic  equator.  Because  of  the  low  air 

3 

density  at  300  km,  the  range  of  keV  electrons  is  very  large  (=»  2.7  xio  km), 
and  it  is  necessary  to  go  to  lower  electron  energies  in  order  to  have  a 
higher  volume  deposition  rate.  The  results  presented  below  are  for  a 325  eV 
electron  beam. 


The  code  results  show  that  the  deposition  volume  can  be  described 
approximately  as  a cylinder  320  km  long  by  30  m in  diameter,  for  a total 
volume  of  2.25  cm^. 


PERMITTED  TRANSITIONS 

The  energy  deposited  in  the  various  states  of  atomic  oxygen  depends 
on  the  ambient  electron  density.  The  higher  the  ambient  ionization,  the  more 
energy  is  imparted  to  it  by  the  electron  beam,  and  the  less  energy  is  available 


for  atomic  excitation.  Table  7-1  shows  the  computed  emission  in  several 
permitted  lines  of  atomic  oxygen  for  three  different  values  of  the  ambient 
election  density.  The  results  in  the  last  six  columns  are  for  a beam 
current  of  1 ampere,  and  the  intensity  is  in  Rayleighs. 

The  intensities  shown  in  Table  7-1  were  computed  for  a path  length 
of  30  m,  corresponding  to  the  diameter  of  the  cylindrical  deposition  volume. 
This  is  appropriate  for  a viewing  geometry  that  includes  a ground-based 
interferometer  located  near  the  magnetic  equator  where  the  long  axis  of  the 
cylinder  is  parallel  to  the  earth's  surface. 

The  intensities  in  Table  7-1  need  to  be  compared  with  the  airglow 
background,  at  the  corresponding  wavelengths,  that  arises  from  radiative 
recombination  of  0^  in  the  F layer.  For  estimating  the  backgrounds,  we 
have  used  radiative  recombination  values  quoted  by  Tinsley  et  al.  (Reference 
27)  for  an  electron  temperature  of  2000  °K.  These  values  are  shown  in 
Table  7-2. 

The  volume  emission  rate  is  given  by  the  equation 

• 2 -3-1 

<f»  = an^  (photons  cm  sec  ) , (7-1) 

and  the  brightness,  for  a path  length  L,  is  then 

B = Lan^^  (7-2) 

provided  the  medium  is  optically  thin. 

For  a vertical  path  length  of  200  km,  the  computed  airglow  inten- 
sities, for  the  wavelengths  corresponding  to  those  in  Table  7-1,  are  shown 

O 

in  Table  7-3.  Intensities  as  high  as  200  R at  7774A  in  the  equatorial 
anomaly  region  have  been  reported  (Reference  28) , The  values  shown  in  the 
last  column  of  Table  7-3,  are  not  atypical  of  what  the  background  intensity 
is  likely  to  be  in  an  electron  beam  experiment  near  the  equator. 


Photon  emission  rate  in  various  permitted  lines  of  atomic  oxygen  at 
300-lcm  altitude  by  a 325  eV  electron  beam  of  unit  ampere  current. 


s. 

p 

CM 

o 

1 

v 

P" 

00 

oo 

to 

r** 

• 

• 

• 

• 

• 

CM 

CO 

1 0)  ^ lO 

1 Q.Q£ 

£ 

CM 

U 

4-) 

fM 

1 

Q) 

<0  to 

CSJ 

o 

CO 

a. 

CL  O 

• 

P“ 

CsJ 

CO 

to 

<j\ 

1 >>  E 

1 +* 

1 

o 

1 (/>  CO  IT) 

c 

o 

CM 

1 

p“ 

r“ 

p~ 

4^ 

n 

VO 

CO 

C 

• 

• 

• 

• 

OJ 

VO 

CM 

c 

♦ 

* 

to 

o 

CO 

CSJ 

CM 

CO 

p* 

CO 

to 

CM 

o. 

• 

• 

• 

E 

a\ 

^mm 

CM 

to 

r- 

u 

o 

in 

VO 

CM 

CM 

CO 

o 

CO 

o 

00 

p“ 

• 

♦ 

* 

'e 

a 

p^ 

r“ 

r— 

CO 

V) 

to 

e 

o 

o 

p^’ 

CO 

r“ 

CM 

CO 

4-» 

II 

CNJ 

o 

to 

o 

« 

• 

• 

* 

sz 

<u 

00 

CM 

CM 

Q. 

c 

“O 

VO 

CO 

tr> 

to 

a? 

o 

1 

1 

1 

1 

1 

p— 

VO 

a> 

00 

•r“ 

• 

« 

V) 

to 

00 

P-* 

CM 

o 

to 

a. 

a; 

a 

CO 

to 

VO 

VO 

1 

> 

VO 

1 

1 

1 

1 

Oi 

o 

VO 

CM 

r— 

00 

fmm 

CM 

CO 

Q) 

Q. 

(/> 

C 

o 

tn 

o 

CO 

1 

1 

VO 

1 

to 

1 

o 

II 

00 

00 

CO 

c 

• 

« 

o. 

CO 

CM 

CM 

to 

00 

■•J 

c 

VO 

00 

VO 

0) 

o 

to 

VO 

r*^ 

P-  I 

CO 

CO 

00 

r*** 

V 

^ 1 

P“ 

oo 

(0 

1 

3 

1 

VO 

00 

1 

to 

• 

00 

CO 

1 

rs.. 

CO 

00 

to 

p— 

1 

VO 

o 

CO 

to 

VO 

CO 

T 

to 

CO 

CO 

• 

CO 

1 

o 

CO 

CM 

Table  7-3.  Computed  airglow  intensity  for  200-km  vertical  path  and 
selected  values  of  ambient  electron  density. 


Wavelength 

(A) 

ne  = 

10  ^ 

Intensity  (R) 

10^ 

5x10® 

1304 

4.8’^ 

4.8 

1.2^ 

1356 

9.6"^ 

9.6 

2.4^ 

4368 

1.0'^ 

1.0‘^ 

2.5 

1 

7774 

6.8'^ 

6.8 

CM 

8446 

3.0"^ 

3.0 

7.5^ 

I A comparison  between  the  last  column  in  Tables  7-1  and  7-3  shows 

j that  for  an  electron  beam  current  of  1 ampere,  the  airglow  background  at  each 

; wavelength  is  likely  to  be  considerably  higher  than  the  intensity  generated 

by  the  beam.  Only  if  we  could  increase  the  beam  current  to  values  between 

10  and  100  amps  could  the  signal  to  background  ratio  become  acceptable.* 

However,  such  currents  for  an  electron  energy  of  325  eV  appear  to  be 

I unreasonably  large  because  of  space-charge  limitations.  For  example,  in 

Project  EXCEDE,  the  typical  beam  current  attained  at  a 3-keV  energy  is  10  amps. 

Assuming  the  current  is  space-charge  limited,  so  that  the  current  varies  as 
3/2 

(voltage)  , then  for  a potential  of  325  volts  the  current  is  only  about 
0.4  amps. 

The  conclusion  is  that  electron  beam  excitation  of  permitted  lines 
of  atomic  oxygen  at  300-km  altitude  provides  too  weak  a source,  relative  to 
the  night  glow  background,  for  measurement  of  doppler  shifts. 

* Even  this  assumes  that  the  electron-excited  emission  volume  completely 
fills  the  sensor  field  of  view,  which  is  unlikely  because  a 30  m region 
viewed  from  a range  of  300  km  would  require  an  angular  field  of  only 
6 10“3  degrees. 


FORBIDDEN  TRANSITIONS 

We  have  looked  at  one  other  line  of  atomic  oxygen,  namely  the 
13  ® 1 

forbidden  transition  D-  P at  6300A.  The  D state  has  a lifetime  of  110 
sec  and  so  the  emitting  region  may  consist  of  an  afterglow  trail  considerably 
longer  than  the  30-m  diameter  of  the  deposition  cylinder.  Consider,  for 
example,  the  geometrical  arrangement  depicted  in  Figure  7-1.  Here,  the 
rocket  trajectory  lies  in  a plane  normal  to  the  earth's  magnetic  field. 

The  instantaneous  deposition  volume  is  still  a cylinder  along  the  field 
lines  (normal  to  the  page),  but  the  emitting  region  will  be  the  volume  swept 
out  by  the  cylinder  as  it  moves  down  along  the  rocket  trajectory.  The 
effective  path  length  through  the  region,  as  seen  by  a sensor  looking  up  the 
rocket  trajectory,  is  the  distance  the  rocket  moves  during  some  fraction 
(>.0.5)  of  the  radiative  lifetime  of  0(^D)  . 


Earth's  Field 
(Normal  to  Page) 


Rocket  Trajectory 


Instantaneous  Deposition 
Volumes, Forming  Trail 
at  6300A 


Sensor  On  Ground  Looking 
UP  Trajectory 


Figure  7-1,  Depicting  geometrical  arrangement  for  viewing  6300A 
emission  from  ground. 
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We  estimate  the  brightness  seen  from  the  ground  as  follows.  If  P 
is  the  production  rate  of  0{^D)  by  the  electron  beam,  then  in  the  absence 
of  quenching,  which  is  negligible  at  300-km  altitude, 

dLQiD)]  , p . a[0(^D)]  (7-3) 

where  A is  the  radiative  transition  probability  (=*=1/110  sec”^).  Thus, 
while  the  beam  is  on, 

[0(^D)I  = J (l-e'^^^  . (7-4) 

But,  because  of  rocket  motion,  any  volume  element  will  be  irradiated  for  only 
a short  time.  At,  during  which  the  concentration  of  0(^D)  will  build  up  to 
the  maximum  value  PAt  and  then  decay  by  radiation  according  to  the  equation 

[0(^D)]  = PAte"^^  . (7-5) 

As  an  approximation,  assume  that  a volume  element  at  a maximum  distance,  L, 
behind  the  rocket  still  contributes  to  the  emission  when  the  radiation  from 
it  is  down  by  a factor  of  2 from  its  peak;  that  is,  when  e”^^^  = 1/2.  For 
1/A=110  sec,  the  time,  t^,  required  to  satisfy  this  condition  is  76  sec. 

The  effective  path  length  through  the  emitting  region  is 

L = vtj  (7-6) 

where  v is  the  rocket  velocity,  and  so  the  brightness  is 

B =*=  j Av  tj  PAt  . (7-7) 

But  At,  the  irradiation  time,  is  the  time  taken  for  the  rocket  to  traverse 
the  diameter,  D,  of  the  deposition  cylinder.  That  is 

At  = ^ . (7-8) 
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Thus 


Ba>-jADPtj^  (photons  cm"^  sec"^)  . (7-9) 

For  A =1/110,  D = 3xio^,  = 76,  we  get 

B « 1 X 10"^  P (Rayleighs)  . (7-10) 

Eqjjation  7-10  was  used,  along  with  the  production  rates  of  0(^D) 
computed  for  three  different  values  of  the  ambient  electron  density,  to 

O 

calculate  the  brightness  of  6300A  emission.  The  results  are  shown  in 
Table  7-4. 

The  intensities  shown  in  Table  7-4  are  to  be  compared  with  the 
tropical  night  glow  at  6300A,  typical  values  for  which  are  400  to  800 
Rayleighs  (Reference  27) . It  is  clear,  therefore,  that  the  beam-induced 
emission  is  small  compared  with  the  night  airglow. 

Table  7-4.  Estimated  intensity  of  6300-A  emission  excited  by  a 325  eV 
electron  beam  of  unit  ampere  current  at  300-km  altitude. 


Wavelength 

(A) 

1 

1 I 

Production  Rate  of  0(^D) 

-3  1 -1 

(cm  sec  amp  ) j 

5 6 6 

n„  = 10'’  10°  5x10° 

e 

Intensity  Per  Ampere 
(Rayleighs) 

n„  = 10®  10®  5x10® 

e 

6300 

B 

^ 1.2^  4.2^ 

7.3^  1 2^  4.2 

Conclusions 

We  conclude  that  the  signal  (from  both  permitted  and  forbidden 
transitions  in  0)  attainable  in  an  electron  beam  experiment  at  300-km 
altitude  is  too  small,  relative  to  the  airglow  background,  to  provide  a 
usable  means  for  measuring  the  neutral  wind  conqjonent  normal  to  the 
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APPENDIX  A 

PARTITION  OF  ELECTRON  ENERGY  (>  1 keV)  IN  ATMOSPHERIC  GASES 


PRELIMINARY 

In  the  discussion  o£  electron  temperatures,  given  in  Section  4 
of  this  report,  or  in  the  calculation  of  the  effects  of  X-ray  or  electron 
deposition,  one  needs  a simple  semi -quantitative  method  of  partitioning  the 
deposited  energy.  In  this  appendix  we  use  the  continuous  slowing  down 
approximation  (Reference  3)  to  develop  approximate  relationships  for  the 
partition. 

PARTITION  FOR  SMALL  FRACTIONAL  IONIZATION 

The  equation  that  determines  the  slowing  down  of  an  electron  in  a 
gaseous  mixture  of  N^,  ©2,  0,  N2*,  ...  is  given  by 

f - -Wzll-N^CE)  - lO^lLg^CE)  - [N)L^(E) 

-[0]Lq(E)  - [N2^]Ljj^^(E)  ...  (A-1) 

where  [N2] , [O2] , . . . , refers  to  the  concentrations  of  N2,  O2,  . . • , and  L refers 
to  their  respective  loss  functions  (References  1,  3).  The  loss  function  is 
defined  by  the  equation 


) 


If 


L = 


z 


o.w.  + ionization  terms 
1 } 


(A-2) 


A-1 


where  a^(E)  is  the  cross  section  for  an  electron  of  energy  E to  excite 
state  j of  the  atom  (or  molecule)  and  lose  an  amount  of  energy  in  the 

process,  E is  the  mean  energy  of  the  electron,  and  s is  the  mean  path 
length  traveled. 

It  is  an  interesting  fact  that  the  above  loss  functions  for  keV 
electrons  are  related  approximately  by 


L»  * 2L  « 2L-  * ^ - 2L 


(A-3) 


That  is,  the  loss  function  depends,  approximately,  only  on  the  number  of 
atoms  present  in  the  gas.  Thus,  we  can  write  Equation  A-1  as 


--AL(E) 


(A-4) 


where  A is  defined  by  the  relation 


A = 2[N2]  + 2[02]  + [0]  + [N]  + ... 


(A-5) 


For  keV  electrons,  L is  given  by 


L(E)  = (keV  cm^) 


Now  the  probability  that  a single  electron  (unit  flux)  of  energy 
E,  traversing  a distance  ds,  will  excite  the  state  of  is  given  by 


INjlCj  j ds 

^ j 


(A-7) 


Therefore,  the  number  of  excited  states  formed  by  an  electron,  per  unit 
distance,  while  it  is  slowing  down  is 


^ j 


(A-8) 


Integrating  Equation  A-8  over  s,  and  making  use  of  Equation  A-4,  we  find 


(A-9) 


Since  ^ 2L,  Equation  A-9  can  be  rewritten  as 


2[N2l 


II 


o 2^ 


(A-10) 


where  E^  is  the  initial  electron  energy  and  the  term  in  brackets  depends 
only  on  the  results  calculated  (or  measured)  in  a pure  gas. 

For  keV  electrons, 

? °[N,1 


(A-11) 


where  e, 


where 


is  a constant. 


Equation  A-10  can  now  be  written  as 


'^[N2]  " ^[N2]''[N2] 

j j 


ffN^]  ^ 2(N2l/A 


CA-12) 


(A-13) 


The  total  amount  of  energy  going  into  state  [N2]j  is  J[N2]. 
'[^2].*  fraction  of  the  total  energy,  E , going  into  [N2].-  is^ 

'’iM  1 “ I ^ 1 '"'rw  1 • (A-14) 


" ^[N2]^[N2]  "'[N2] 

j j 


A,.  ■ I ^ j|-|  - , - 


Clearly,  ^[N2]-”[N2]-  partition  one  would  calculate  for  a pure 

gas.  These  results  can  be  applied  to  any  gas  in  the  mixture.  Equation  A- 14 
tells  us  that  we  need  only  combine  the  partition  calculated  in  pure  gases 
linearly,  weighting  each  by  the  fractional  number  of  atoms  present. 

The  above  derivation  was  for  the  primary  electrons  only,  but  when 
secondaries  and  higher  order  effects  are  included,  detailed  calculations  show 
that  Equations  A-10  through  A-14  are  still  approximately  correct  although 
the  constants  e are  then  redefined.  The  fractional  amount  of  energy  going 
into  each  excited  state  is  approximately  independent  of  the  initial  electron 
energy.  In  particular,  this  holds  for  states  of  ionization.  For  example, 
it  is  well  known  that  the  number  of  ion-pairs  formed  per  eV  deposited  in 
"normal  air"  is  a constant  over  a wide  range  of  incident  electron  energies. 

In  the  present  notation  this  can  be  written  as 


ion-pairs/eV  « Y]  f.e. . 

" 1 ij 

11 


(A-15) 


where  i denotes  the  gas  and  j denotes  a state  of  ionization. 

In  Tables  A-1  to  A-4  are  listed  the  constants  e for  N^,  0^, 

N,  and  0 gases.  These  results  include  all  orders  of  secondary  electrons. 
The  actual  solution  is  an  integral  equation  which  can  be  solved  exactly 
(Reference  15).  Since  the  above  derivation  is  approximate  rather  than  exact, 
we  have  weighted  the  results  so  that  N2  and  O2  reproduce  the  results  of 
an  80  percent  N2  percent  0^  gas  mixture  while  0 is  weighted  to 

a 50  percent  N2  and  50  percent  0 gaseous  mixture. 

These  results  can  now  be  used  to  calculate  ionization,  excitation 
and  partition  of  energy  in  a mixture  of  gases  following  X-ray  or  electron 
deposition.  The  effect  of  the  variation  of  atmospheric  constituents  during 
deposition  can  also  be  taken  into  account  in  a straightforward  manner.  For 


Table  A-1.  Partition  In  a pure  N2  gas. 


i 

State 

^[Ngl 

Mean  Excitation  Energy 
(eV)* 

Energy  Partition 
(Pi) 

1 

.0252 

16.98 

.427 

2 

n‘^(^p) 

.002 

19.42 

.039 

3 

N^(^D) 

.0008 

21.32 

.017 

4 

N^(^S) 

.0002 

23.52 

.005 

5 

N(^S) 

.0078 

4.88 

.038 

6 

N(^D) 

.007 

7.26 

.051 

7 

N(^P) 

.0018 

8.46 

.015 

8 

N2(v)direct 

.223 

.3 

.067 

9 

N2(v)cascade 

.007 

.3 

.002 

Ion  Kinetic  Energy 

.008 

Neutral  Kinetic  Energy 

.023 

Energy  Into  Electron  Gas  and  N2  Rotational  Excitation  .031 

Energy  Into  Radiation  and  Metastable  States  of  N2 

.277 

* The  mean  energy  is  not  necessarily  the  energy  of  the  lowest  lying 
state  since  several  states  may  have  been  combined  to  form  the 
above  table.  For  atomic  products,  one  half  the  dissociation 
energy  has  been  assigned  to  the  mean  excitation  energy. 
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Table  A-2.  Partition  In  a pure  O2  gas. 


1 

State 

*"[02] 

Mean  Excitation  Energy 
(eV) 

Energy  Partition 
(P^) 

10 

O2' 

.0203 

16.75 

.335 

11 

O^S) 

.007 

16.15 

.113 

12 

0‘^(^D) 

.004 

19.45 

.078 

13 

0(^P) 

.0387 

2.54 

.098 

14 

1 

0(^0) 

.0127 

4.5 

.057 

15 

02(v) 

.67 

.25 

.168 

16 

02(^A) 

.040 

Ion  Kinetic  Energy 

.002 

Neutral  Kinetic  Energy 

.031 

Energy  Into  Electron  Gas  and  O2  Rotational  Excitation  .005 

Energy  1>nto  Radiation  and  Metastable  States  of  Og 


.073 


Table  A-3.  Partition  In  a pure  0 gas. 


r'  — 

1 

State 

^[0] 

Mean  Excitation  Energy 
(eV) 

Energy  Partition 
(Pi) 

17 

0(^0) 

.110 

1.96 

.217 

18 

0(^S) 

.005 

4.98 

.021 

19 

O^(^S) 

.016 

13.61 

.218 

20 

O'^(^D) 

.014 

16.9 

.237 

21 

0‘^(^P) 

.007 

18.5 

.130 

0 

Energy  Lost  to  Electron  Gas  and  Fine  Structure  of  0(  P)  .058 

Energy  Into  Radiation 

.119 

Table  A-4.  Partition  In  a pure  N gas. 

1 

State 

^[N] 

Mean  Excitation  Energy 
(eV) 

Energy  Partition 
(Pi) 

22 

N(^D) 

.082 

2.38 

.195 

23 

N(^P) 

.022 

3.58 

.079 

24 

N^(^P) 

.018 

14.5 

.261 

25 

-016 

16.4 

.262 

26 

.007 

18.6 

.130 

Energy  Lost  to  Electron  Gas  .008 


Energy  Into  Radiation 


.065 


example,  if  S represents  the  total  amount  of  energy  deposited,  the  rate 
at  which  [N2]  is  destroyed  is 


'‘W2I  ,1  rr 

dt  - - A jh  ■ 


(A- 16) 


where  the  values  of  are  given  in  Table  A-1, 


In  terms  of  the  variable  U = «?/A,  Equation  A-16  can  now  be  written 


‘i[N2l 


■ (A-17) 


with  solution 


where 


(N2]  = [N2]  e 

0 


R,  = 2e,  + e-  + e_  + e.  + + e,  + e_ 

1 1234567 


(A-18) 


and  [N_]  is  the  initial  concentration  of  N_  at  the  start  of  the  depo- 
z o z 

sition.  Concentrations  of  other  states  can  be  calculated  using  a similar 
procedure . 

PARTITION  FOR  LARGER  FRACTIONAL  IONIZATION 

Results  such  as  A-18  are  valid  provided  the  fractional  ionization 
remains  small.  If  it  does  not,  then  we  must  add  an  additional  term  to  the 
right  side  of  Equation  A-4,  rewriting  it  as 


- = - AL  - n L 
s e e 


(A- 19) 


1 


I 


;i 


> ( 


i 


Here,  n is  the  electron  density  and  L is  the  electron  loss  function 

C 6 

given,  for  keV  electrons,  by  (Reference  15) 


Lg  =<  3.4  X 10‘^Ve  keV  cm^ 


(A- 20) 


function. 

get 


In  general,  the  partition  of  energy  now  becomes  a very  complex  ^ 

Factoring  A out  of  the  right  hand  side  of  Equation  A-19,  we  | 


g . - A(L  . (A-21) 

where 

fg  = n^/A  . {A-22) 


The  fraction  of  the  energy,  E^,  going  into  the  electron  gas  is  given  by 


P 

e 


%(Ec 


•I) 
n(E,T)n 


L (T) 
e e^ 


AL(T)  *n^Lg(T) 


dT 


(A-23) 


o o 

where  n(E,T)  is  the  number  of  secondary  electrons  of  energy  T produced 
by  a primary  of  energy  E,  and  I is  the  ionization  potential  (see 
Reference  15) . 


We  see  that  as  long  as  f L is  small  compared  to  L,  the  previous 

© c 

results  hold..  When  f L is  large  compared  to  L,  then  all  of  the  energy 

c © 

will  go  into  the  electron  gas.  The  region  in  between  is  quite  complex 

particularly  since  the  partition  depends  on  an  integral  over  energy,  and  L 

and  L have  different  energy  dependencies.  Thus,  primary  and  secondary 
© 

electrons  behave  differently  as  a function  of  fractional  ionization,  f . 

© 

To  first  order,  we  can  consider  four  separate  regions  depending  on  whether 
it  is  primary  or  secondary  electrons  that  are  mainly  involved  in  heating  the 
electron  gas. 

A-9 


I 


We  first  note  that  typical  energies  of  secondary  electrons  are 
-10  eV  while  primary  electrons  are  100  eV  to  10  keV  or  more.  Thus,  L 
for  a secondary  electron  is  several  orders  of  magnitude  larger  than  L 

Q ® 

for  a primary.  For  "small"  fractional  ionization  (fg<10  ) the  principle 
method  of  heating  the  ambient  electron  gas  is  by  secondary  electrons  and 
essentially  no  heating  is  done  by  the  primaries.  The  heating  rate  is  so 
small  that  it  can  be  considered  as  essentially  a perturbation  on  the  previous 
results.  Thus,  in  Region  1 the  first  term  on  the  right  hand  side  of  Equa- 
tion A-23  is  zero,  and  the  partition  of  energy  into  the  electron  gas  becomes 

’sCEo-l) 


1 f n(E.T)n^L^ 
“ J AL(T) 


dT  = f 6, 
e 1 


(A- 24) 


where  6^  is  a constant  applicable  to  Region  1. 

As  the  fractional  ionization  increases,  however,  there  comes  a 
point  where  f L > L for  secondaries  but  f L « L for  primaries . That 
is,  we  reach  a point  of  saturation  where  all  of  the  secondary  electron 
energy,  but  none  of  the  primary  energy,  goes  into  heating  the  ambient  electron 
gas.  In  this  region  the  partition  of  energy  into  the  electron  gas  becomes 
simply  a constant,  independent  of  the  fractional  ionization.  Thus,  in 
Region  2, 


Pe«62 


(A- 25) 


Region  3 is  where  f L <L  for  primaries,  but  the  primaries  still 

6 6 

contribute  to  heating  the  ambient  electron  gas.  The  first  term  of  Equation 
A-23  can  then  be  approximated  as 


‘'o 

2.  r V;! 
=oj 


dE  * f 


(A- 26) 


The  partition  then  becomes,  approximately. 


’’e  " " ^e  ^3 


(A-27) 


where  the  term  £^62^^  comes  from  the  reduced  amount  of  energy  going  into 
secondary  electrons. 


Finally,  Region  4 is  where  f L »L  for  all  energies 


and  P 1 . 
e 


The  foregoing  behavior  is  illustrated  schematically  in  Figure 
A-1.  We  see  that  in  Regions  1 and  3,  varies  linearly  with  f while 
in  Regions  2 and  4 it  is  constant. 


Region  3 Region  4 


Figure  A-1.  Schematic  illustration  of  the  dependence  of 

P-  on  r . 
e e 


Detailed  runs  of  the  ARCTIC  code  (References  1,  3)  tend  to 

it 

confirm  this  analysis.  We  have  done  a parametric  study  over  a wide  range 
of  energies  (1<E^<1S  keV)  and  over  compositions  corresponding  to  a variety 
of  atmospheric  mixes  between  60  and  200-km  altitude. 

Figure  A-2  presents  our  results  for  one  energy  (E^  = 4 keV)  at  two 
compositions  corresponding  to  altitudes  of  100  and  184  km.  We  see  that  the 
results  are  quite  independent  of  the  incident  electron  energy.  Regions  2 
and  3,  however,  tend  to  overlap,  giving  rise  to  an  intermediate  behavior. 

The  partition  for  fractional  ionizations  f > lO"^  varies  approximately 

25  ® 

as  f ■ .As  expected,  it  is  nearly  independent  of  atmospheric  compo- 

® _g 

sition.  In  Region  1 (f  < 10  ) , the  variation  of  P with  f is  indeed 

seen  to  be  linear,  but  the  partition  is  composition  dependent.  This  is 
not  surprising  in  view  of  the  low  energies  of  the  secondaries  (~  10  eV) . 


A-12 


Altitude  = 184  km 


Initial  Fractional  Ionization  (fg) 

Calculated  fractional  energy  deposited  in  electron  gas  as  function  of 
initial  fractional  ionization. 


Y2  = 1-75  X 10^ 

M = [Molecules] /([Molecules]  + 12  x [ATOMS]) 


To  first  order,  the  complete  partition  of  energy  into  all  modes 
is  given  by  Equation  A-28  (for  the  ambient  electron  gas)  and  by  the  results 
of  Tables  A-1  through  A-4  reduced  by  a factor  (1  - Pg) • 

The  modifications  presented  above  will  probably  not  affect  calcu- 
lations below  altitudes  of  100  km.  However,  above  200  km,  where  the 
fractional  ionization  regularly  exceeds  10”°,  the  changes  can  be  dramatic. 
In  particular,  use  of  Equation  A-24,  rather  than  Equation  A-28  at  these 
altitudes  (the  assumption  usually  made),  will  produce  erroneous  results. 


APPENDIX  B 

PARTITION  OF  ELECTRON  ENERGY  (<  25  eV)  IN  ATMOSPHERIC  GASES 


The  partition  of  low  energy  electrons  (<  25  eV)  among  various 
states  is  often  needed  for  chemistry  or  IR  calculations.  These  low  energy 
electrons  can  arise  from  UV  deposition  or  from  deexcitation  of  excited  atomic 
and  molecular  states.  We  have  performed  a parametric  study,  using  the  ARCTIC 
code,  of  this  partition  as  functions  of  electron  energy  and  fractional  ioni- 
zation. The  results  are  presented  here. 

Figures  B-1  through  B-4  show,  for  deposition  in  N^,  0^,  N,  and 
0,  respectively,  the  fraction  of  the  electron's  energy  deposited  in  the 
thermal  electron  gas  + rotational  states  + fine  structure  states  of  0(^P) 
as  functions  of  the  initial  fractional  ionization,  f^,  (see  Equation  A-22) • 

We  make  the  assumption  that  initially  all  of  this  energy  is  deposited  in 
the  final  (after  deposition)  thermal  electron  gas.  This  is  because  the 
energy  going  into  rotation  or  into  fine  structure  excitation  comes  about 
only  after  most  other  excitations  have  occurred,  at  which  time  the  electron 
has  lost  most  of  its  energy.  The  cross  sections  for  these  processes  are 
very  flat  at  low  energies  so  that  subsequent  cooling  of  the  electron  gas  by 
rotational  or  fine  structure  excitation  will  be  approximately  correct. 

T 

; Over  the  energy  range  considered  (0-24  eV)  the  form  of  all  the 

I curves  shown  in  Figures  B-1  through  B-4  can  be  represented  approximately  by 

\ Pg  = (A/E)“  (B-1) 

B-1 

I 

■i 

A 


Electron  Energy  (eV) 

B-1.  faction  of;;, «ro"f  r^Tr.Sl^'loint’tfonTn  Ture  Nj 


Fracti on 


where  P is  the  fraction  of  the  hot  electron's  energy  going  into  the 
electron  gas,  E is  the  initial  energy  of  the  hot  electron,  and  A and  a 
are  parameters  that  depend  on  the  fractional  ionization.  These  parameters 
are  given  approximately  by 


a »MAX[0.14,  MIN(1,  0.96-8x10^  f^)] 


(B-2) 


where 


A « 24“(U  + EXP[log(l-U)  (14.5/(14.5  + log  f^))^'^^])  (B-3) 


U = MAX[.075,  8.616x  10"  T^/E] 


In  addition,  we  have  calculated  the  energy  going  into  other  states. 
Figures  B-5  and  B-6  show,  respectively,  the  fractions  of  energy  going  into 

vibrational  states  of  N.  and  0 for  fractional  ionizations  less  than 

•6  ^ * 
lO"  . We  have  approximated  these  curves  by  the  functions 


N2(v) 


^ (1  - 2/E)  E > 2 


(B-4) 


(1  - 1/E)  E > 1 


02(v) 


(B-5) 


These  approximations  are  plotted  in  Figures  B-5  and  B-6. 

In  addition,  we  have  calculated  the  number  of  atomic  nuclei  formed 
per  eV  deposited,  as  was  done  in  Appendix  A.  These  results,  shown  in 
Figures  B-7  through  B-10,  have  been  approximated  analytically.  The  results 
are  shown  in  Table  B-1. 


Electron  Energy  (eV) 


Figure  B-6 


. Fraction  of  electron  energy  deposited  In  vibrational 
excitation  of  0«  for  zero  fractional  Ionization. 
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Electron  Energy  (eV) 

Figure  B-10.  Number  of  0(^0)  atoms  formed  per  electron  deposited  in  1 

0 gas.  I 

I 
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Table  B-1.  Species  formed  in  low  energy  electron  deposition. 


Gas 

Atomic  Product  Formed 

Number  Created  Per  eV  Deposited 

^2 

N(S) 

0.39(1-10/E)^'®  ; 10<  E<  24  eV 

N(^0) 

0.38(1-14/E)^*^  ; 14<E<24  eV 

N(^P) 

0.17(1-14/E)^-^  ; 14<  E^  24  eV 

O2 

0(^P) 

0.09(E-5.4)  ;5.4<E<24eV 

0(^0) 

0.07(E-9.5)  :9.5<E<24eV 

N(^S) 

N(^D) 

0.6  ; 4<E<24eV 

N(^P) 

0.22  ; 6<E<24eV 

0(^P) 

0(^0) 

9.8E"^1-1.96E"^)^;1.96<E<24  eV 

For  mixtures  of  gases,  we  have  assumed  linear  superposition  using 
the  weighting  factors  of  Appendix  A (Equation  A-12) . In  addition,  when  the 
fractional  ionization  exceeds  10~^,  the  fractions  shown  in  Figures  B-5 
through  B-10  must  be  reduced  to  allow  for  the  additional  energy  going  into 
ionization.  For  this  purpose,  we  have  used  a reduction  factor  g^,  where 


g^  = (l-y/(l-U) 


(B-6) 


Thus,  the  partition  of  energy  into  and  O2  vibration  can  now  be 

written  as 


and 


'’n2(v) 

’’02(v) 
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E > 2 eV 


E > 1 eV 


(B-7) 

(B-8) 
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From  Figures  B-1  through  B-10  we  see  that  our  formulas  are  approxi- 
mate rather  than  exact.  Therefore,  it  may  occur  that  energy  is  not  strictly 
conserved.  We  have  corrected  for  this  by  multiplying  the  results  by  a 
normalization  factor,  R,  defined  by 


R = E/(sum  of  all  energy  lost) 


(B-9) 


APPENDIX  C 

ON  OPTICAL  MEASUREMENT  OF  AURORAL  NO  PROFILES 


INTRODUCTION 

We  have  performed  a study  to  determine  the  feasibility  of  measuring 
the  enhanced  NO  profile  in  an  auroral  environment  by  rocket-borne  photometer 
observations  of  the  continuum  emission  from  the  N0  + 0-*-N02  + hv  reaction. 

The  work  was  originally  motivated  by  a requirement  for  additional  experiments 
to  be  carried  on  board  a multi-instrumented  rocket  flown  into  the  aurora. 


The  idea  is  basically  to  measure,  as  a function  of  altitude,  the 
continuum  emission,  in  a specified  wavelength  band,  with  a horizontal-viewing 
photometer.  These  measurements,  combined  with  in  situ  data  on  the  0 con- 
centration (from  a mass  spectrometer)  and  the  known  2-body  rate  constant  for 
the  continuum  reaction,  should  be  sufficient,  under  certain  conditions,  to 
infer  an  NO  profile. 

The  results  of  the  study  are  presented  in  this  appendix.  They 
provide  two  specific  wavelength  bands  in  which,  for  an  enhanced  NO  con- 
centration, the  continuum  emission  should  sufficiently  dominate  the  back- 
ground noise  to  render  the  experiment  feasible. 

The  basic  problem  is  to  find  suitable  wavelength  regions  for  which 
the  continuum  emission  will  not  be  masked  by  background  noise  from  starlight  + 
zodiacal  light + diffuse  galactic  light + auroral -induced  air  fluorescence. 

This  problem  is  first  addressed. 


CONTAMINATION  RADIATION 


Air  Fluorescence 

Unfortunately,  most  wavelength  regions  of  interest  for  this  experi- 
ment are  profusely  populated  by  air  lines  and  bands  that  can  be  excited  by 
the  drizzle  precipitation.  Two  regions  that  are  relatively  clean,  however, 

o o o o 

are  (5025A-5075A)  and  (4750A-4850A) . This  can  be  seen,  for  example,  from 
References  29  and  30,  and  from  a survey  of  atomic  and  molecular  band  tables. 
It  is  for  these  two  spectral  regions  that  we  have  done  our  calculations. 

The  following  table  shows  the  principal  bands  that  fall  within  the  two 


regions . 


Table  C-1.  Fluorescence  bands  in  two  wavelength  regions. 


Wavelength  Region 


Interfering  Bands 


5025A-5075A 


N,(1P;  10,  4;  5068A) 
N^CIP;  11,  5;  5045A) 
N2(1P;  12,  6;  5027.A) 
N2(2P;  0,  6;  5032A) 
N2(2P;  3,  10;  5067A) 
N2(VK;  7,  19;  5030A) 


N^CVK;  0,14;  5061A) 

N2^(1N;  2,  5;  5077A) (degraded  to  violet) 


4750A-4850A 


N2(VK;  2,  15;  4837A) 
N2(VK;  5,  17;  4770A) 
N2(2P;  2,  8;  4813A) 
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We  have  calculated  the  volume  emission  rate  from  these  bands  using 
the  following  formulae: 


*(lP;v"v'")  = 


2,9  q': 
^v  o 


q'  q(.^X 
VO 


3-r  nN2(C)]  (C-1) 


$(2P;v"v"')  = 


Q ^ Q > > 
VO  V V 


4(VK;v;v'0 


, 0 

• « ^ LJ  ^ A ^ 

V 


P[N2(C)] 


I.IS  q' 


(C-2) 


--  (1)  ^ 


$ClN;v^v'')  = 


<,-'o 

(>;••:■■)■£  KW) 


(C-4) 


Here  the  q are  Franck  Condon  factors  referring  to  the  following  electronic 
transitions  in  N2  and 

q'(c\  ->  X^Zg)  ; ■>  b\) 

q"CB%-xVg)  ; q^^\xX^xhl) 

A^Z*)  q'"CN2*(B)  ■*•  N2(X)) 


q^^^N2^(B)  N2*(X)) 


and  P[N_(C)]  and  P[N»  (X)]  are  the  volume  production  rates  of  N»(C'*n  ) 

2 + ^ ^ u 

and  N2  (X  Eg  ),  respectively,  by  electron  bombardment.  In  Equation  C-1, 

the  second  term  in  brackets  is  the  contribution  to  IP  radiation  from  C-’-B 

cascade.  In  Equation  C-3,  the  second  term  in  brackets  is  the  contribution 

to  VK  emission  from  B->A  cascade.  Furthermore,  Equations  C-1  and  C-3 

assume  (from  the  ARCTIC  code)  that 

P(N2(A^E)  * 1.15  P[N2(B%)] 

and  _ . 

P [N.,  (B-^n)  ] ^ 2 . 9 P [N2  (C  n) ) 


The  volume  emission  rate  from  the  Vegard  Kaplan  bands,  represented 
by  Equation  C-3,  should  be  reduced  to  allow  for  quenching  of  the  A state 
by  atomic  and  molecular  oxygen.  Specifically,  if  the  quenching  reactions 
assumed  are 

N2(A)^.  + 0 iiili.  N2(X)  + 0 (C-5) 

N^CA)^.  ^2  ^^^^N2(X)  + O2  , (C-6) 


then  the  VK  emission  should  rereduced  by  the  factor,  F,  where 


F = 


( ki(v')[0]  +k2(v')[02]  1'^ 


(C-7) 


Here,  A^^^  is  the  Einstein  transition  probability.  Since  rate  constants 
for  quenching  of  individual  vibrational  states  of  N2(A)  are  not  known, 
we  use  instead  the  overall  rate  constants  kj  and  k2  and  then  approxi- 
mate Equations  C-7  by 


F * 


i 


1 + (k,[0] 


+ k (0  ])T 


1-1 


CC-8) 
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where  t is  the  radiative  lifetime  of  the  A state  («  2 sec) . Values  for 
kj  and  are  taken  equal  to  1.5  x 10'^^  and  6.5x10*^^  cm^  sec”^, 
respectively  (References  il,  32,  2 (footnote  on  p.  108)). 


I 
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The  volume  emission  rate  due  to  air  fluorescence  was  calculated 
using  the  foregoing  equations  along  with  the  volume  production  rates 
P(N2(C)]  and  P[N2*(B)]  determined  by  the  ARCTIC  code  for  the  drizzle 
flux  shown  in  Figure  2-8  of  Reference  2.  The  intensity  of  the  air  fluor- 
escence was  then  calculated  by  integrating  the  volume  emission  rate  along 
lines  of  sight  (spherical  earth),  assuming  an  unlimited  extent  of  the 
drizzle  region. 

Star! ight+ Zodiacal  Light + Diffuse  Galactic  Light 

A summary  of  the  contribution  from  these  sources  is  given  in 

Q 

Reference  33, p.  134.  We  have  used  the  value  0.7  R/A.  This  includes  0.35 

O 

R/A  for  the  mean  sky  background  due  to  faint  stars.  However,  on  March  21 
near  midnight  at  Poker  Flat,  the  galactic  pole  is  at  a zenith  angle  of 
about  52*  toward  the  north,  and  the  galactic  equator  (Milky  Way)  is  below 
the  horizon.  Thus  the  value  of  0.35  R/X  is  probably  too  large,  especially 
for  sight  paths  to  the  north,  and  our  value  of  0.7  R/A  is  likely  an  upper 
limit. 

EMISSION  FROM  NO  + 0 CONTINUUM 


i 


At  altitudes  above  about  95  to  100  km,  the  reaction  between  NO 

and  0 that  produces  the  air  afterglow  continuum  is  essentially  bimolecular 

(References  34,  35)  with  a rate  constant  over  the  bands  of  interest  of  about 
-21  3 -1  ®-l 

1 X 10  cm  sec  A . With  this  rate  constant,  the  volume  emission  rate 
from  the  continuum  was  computed  for  the  two  wavelength  bands  using  the  two 
different  NO  profiles  shown  in  Figure  2-11  of  Reference  2.  One  of  them 
corresponds  to  a "typical”  low  latitude  profile;  the  other,  to  the  profile 
inferred  for  the  1973  Multi  event. 
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The  results  are  shown  in  Figures  C-1  through  C-3.  Figure  C-1 
compares  the  continuum  radiance  along  a horizontal  path  with  the  total 

Q O 

background  due  to  air  fluorescence starlight  etc.  for  the  5025A-5075A 
band.  About  half  of  the  background  emission  is  seen  to  come  from  the  extra- 
terrestrial sources.  The  important  point  to  notice,  however,  is  that  if  the 
NO  is  enhanced  as  much  as  assumed,  the  afterglow  continuum  should  dominate 
the  background  by  an  order  of  magnitude  below  100  km  altitude  and  fall  off 
with  increasing  altitude  much  more  rapidly  than  does  the  background  intensity. 
These  two  effects  should  be  observable.  However,  if  the  NO  is  not 
enhanced,  but  corresponds  more  to  the  low-latitude  profile,  the  continuum 
emission  should  be  masked  by  the  background  which  should  remain  fairly 
constant  in  intensity  with  increasing  altitude  to  at  least  115  or  120  km. 

Figure  C-2  is  a corresponding  comparison,  also  for  the  5025A-5075A 
band,  for  a vertical  sight  path  rather  than  a horizontal  one.  It  is  clear 
from  this  figure  that  the  background  emission  (mainly  extraterrestrial) 
dominates  the  afterglow  continuum  a':  almost  all  altitudes,  even  for  the 
enhanced  NO  profile.  Thus,  zenith  measurements  are  clearly  inappropriate 
for  purposes  of  this  experiment. 

Figure  C-3  is  another  comparison  for  horizontal  sight  paths  but 

e o 

for  the  wavelength  band  4750A-4850A.  Also  shown  in  this  figure  is  the 

O * 

computed  3914A  intensity  along  the  same  horizontal  paths  (this  curve  applies, 
of  course,  also  to  Figure  C-1). 
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Comparison  between  Figures  C-1  and  C-3  shows  that  there  is  no  great 
advantage  of  one  wavelength  band  over  the  other.  The  wider  band  provides 
a greater  intensity  (factor  of  2)  for  the  continuum,  but  the  afterglow 
intensity  drops  to  background  at  a slightly  lower  altitude. 
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IO  + 0-!^N02  + hv 

;(5025A<  A<  5075A)  * 5x10"^°  cm^  sec"^ 


A 


Total  Background 


Starlight  Zodiacal 
+ Diffuse  Galactic 


(d  • N^dP-.io.a)  + NgdPdi.s)  ♦ NgdPdz.e) 

♦ N2(2P;0,6)  ♦ N2{ZP;3,10)  + N2(VK;7,19) 
+ I«2(VK;0,14)  ♦ N2*(1N;2.5) 
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Figure  C-1.  Calculated  intensity  in  5025&S  X i5075A  band  from  N0  + 0-*-N02  + hv 
continuum  compared  with  auroral  drizzle- Induced  fluorescence  from 
N2(VK) +N2(1P) +N2{2P)  bands  + starlight  + zodiacal  light galactic 
light  (horizontarpath). 
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Bands  Included  = NgllPllO.A)  + NgdP-.ll.B)  + NgClPdZ.S) 

+ N2(2P;0,6)  + N^CZPd.lO)  + N2{VK;7,19) 


+ N,(VK;0,14)  ♦ N,  (1N;2,5) 


Altitude  (km) 

o o 

Figure  C-2.  Calculated  intensity  in  5025A  < A < 5075A  band  from  NO+0  -►N02  + hv 
continuum  compared  with  auroral  drizzle-induced  fluorescence  from 
N2(VK) +N2(1P) +N2{2P)  bands  + starl ight  + zodiacal  light  + diffuse 
galactic  light  (vertical  path). 


CONCLUSIONS 

The  foregoing  considerations  suggest  that  the  experiment  has 
merit.  However,  success  would  depend  on  the  availability  of  filters  in  the 
appropriate  wavelength  regions  with  adequate  out-of-band  rejection.  These 
are  items  that  we  have  not  addressed. 
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